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SUMMARY 


Possible contamination of the upper atmosphere from the by-products 
of an industrial society has created the need for regular sampling of 
high-altitude atmospheric components. Atmospheric sampling has been 
carried out by NASA for a number of years using U2 aircraft. These air- 
craft have insufficient flight altitude capability for monitoring the 
growth of some potential contaminants which may be generated by aerosol 
container usage. This report examines the increase in sampling altitude 
which could be obtained if the U2 flights were supplemented by flights 
using an available high-performance supersonic aircraft, the Phantom F4-C. 

Altitude potential of an off-the-shelf F4-G aircraft is examined in 
this report. It is shown that the standard F4-C has a maximum altitude 
capability in the region from 85000 to 95000 ft, depending on the minimum 
dynamic pressures deemed acceptable for adequate flight control. By 
using engine overspeed capability and by making use of prevailing winds 
in the stratosphere, it is suggested that the maximum altitude achievable 
by an F4-C should be in the vicinity of 95000 ft for routine flight oper- 
ation by NASA personnel. This altitude is well in excess of the minimum 
altitudes which must be achieved for monitoring the possible growth of 
suspected aerosol contaminants. 
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There have been recent suggestions that use of freon powered aerosol 
pressure containers can lead to significant modifications in upper 
atmospheric constituents over an extended period of time. If this happens 
one of the more serious consequences would be degradation in the ultra- 
violet radiation shielding properties of the upper atmosphere. It is 
thought that aerosol container atmospheric degradation can be detected by 
monitoring the relative density of a variety of chemical specie over a 
period of several years. A prime candidate for such a monitoring process 
is the chlorine radical C^O. Figure 1 prepared by NASA's Ames Research 
Center provides a prediction of the anticipated density increase of this 
atmospheric component over the next sixty years. 

Ames Research Center has now been engaged in upper atmosphere sampling 
for several years. Samples are collected by U2 aircraft. This vehicle 
has a cruise altitude capability of about 70000 feet. Predicted increases 
«n Cj^O concentrations at this altitude during the next sixty years. Figure 1, 
?re small. Therefore, if the growth of this upper atmospheric constituent 
is to be effectively monitored an aircraft having considerably higher 
altitude potential is required. Ideally this aircraft would be an existing 
vehicle, easily maintained, and having an altitude potential on the order 
o£ 90000 feet. Several modern military fighter aircraft have the potential 
for this altitude if a zoom climb maneuver is employed. The zoom climb 
maneuver concept is based on a kinetic potential energy exchange in which 
the vehicle is accelerated to high speed and then exchanges velocity for 
altitude. An extreme form of this maneuver is typically employed when 
establishing aircraft time-to-climb records. For example, the F4 
Phantom fighter has an altitude potential in excess of 100000 
feet if the day, location, and vehicle weight are carefully selected. The 
Phantom is a prime candidate for use in upper atmospheric sampling. It 
has 

a) Wide distribution and availability, over 4000 built in 
several variants. 

b) Demonstrated high altitude capability. 

c) In certain versions, notably the RF4-C, an equipment bay 
large enough to contain the sampling equipment. 

This report defines the altitude potential of standard F4-C aircraft 
using zoom climb maneuvers. Maximum altitude sensitivity to various para- 
meters is presented. The parameter variations considered are 

1) Zoom climb initial weight. 

2) Minimum flight dynamic pressure. 

3) Thrust improvements (available from engine overspeed or 
cold day operations) . 

4) Stratospheric wind distributions. 
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Trajectories are optimized by small perturbation techniques using 
the variational calculus procedure described in the following section. 

A nominal control history. Figure 2, is used to generate a nominal 
flight path. This path supplies the altitude history. Figure 3. The 
control history is then perturbed by a given amount represented by the 
shaded area in Figure 2. Of all control perturbations of given magnitude 
that producing the largest altitude gain. Ah, is employed to generate 
a new trajectory and the process is repeated until further altitude gains 
are impossible. 

The aircraft and planetary representations and simplifications for 
this study are summarized in Figure 3. All trajectory computations in 
the present study were performed using the Atmospheric Trajectory 
Optimization Program, ATOP. This program has been developed with both 
NASA and USAF funding over a period of many years. The program, its 
past applications, and its derivatives are discussed in References 1 
through 12. Program details are given in a later section. For those 
interested in results and not methodology the next two sections should 
be passed over. 
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FIGURE 1. EXAMPLE OF ALTITUDE DEPENDENCE - FREON EFFECTS, 1-D MODEL 


(Terminate in 10 years, 30 year life in troposphere) 
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FIGURE 2. 

SMALL PERTURBATION METHOD OF VARIATIONAL CALCULUS 
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FIGUIU2 3. AIRCRAFT AND PLANETARY MODELS 



T=T(M,h,N); N=1 

Wp=WF(M,h,N);N=l 

CL=CL(M,h,ot) 
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THE STEEPEST DESCENT METHOD 


Problem Statement 

Point mass motion is governed by three second order 
differential equations of position together with a first 
order differential equation governing the mass. By suit- 
ably defining additional state variables, it is possible 
to reduce these equations to a set of first order differ- 
ential equations. Point mass motion is, therefore, governed 
by a set of first order differential equations. The form 
of these equations is 




n - 1,2 ........ N 

D w 1,2 M (i) 


That is, there are N state variables whose derivatives 
Xj^(t) are defined by N first order differential equations 
involving the state variables, together with M control vari- 
ables, , and t, the independent variable itself. 

Constraints may be imposed on a set of functions of the 
state variables and time at the end of the trajectory. In 
this case, a set of constraint functions of the form 

l*p| - ° 

P m 1,2 .......P (2) 

can be constructed which the final trajectory must satisfy. 

Any one of the constraints may be used as a cut-off function 
which, when satisfied, will terminate a particular trajectory. 
The cut-off function can, therefore, be written in the form 



and determines the trajectory termination time T. In all, 
then, when the cut-off function is included, there are 
(P + 1) end constraints. 

Finally, it may be that some other function of the state 
variables and time at the end of the trajectory is to be 
optimized. Hence, a pay-off function 



which is to be maximized or minimized, can be constructed. 
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Now, suppose that a nominal trajectory is available. The 
requirements of this trajectory are modest; it must satisfy 
the cut-off condition, Equation (3) , but it need not optimize 
the pay-off function or satisfy the constraint equations. To 
generate this nominal trajectory by integrating Equations (1) , 
the vehicle characteristics, the initial state variable values, 
and a nominal control variable history must be known. Once 
this nominal trajectory is available, the steepest descent 
process can be applied. To do this, the trajectory showing the 
greatest improvement in the pay-off function, while at the 
same time eliminating a given amount of the end point errors 
as measured by Equations (2) for a given size of control vari- 
able perturbation, is obtained by application of the Variational 
Calculus. 


Equations (2) provide an end point error measure, for they 
will only be satisfied if the end points have been achieved. 
Therefore, any non-zero represents an end point error which 
must be corrected. A convenient measure of the control variable 
perturbation can be defined by the scalar quantity, 


Dp2 


C L‘“H[ 


w(t) 


ia(t)> dt 


(5) 


where W is any arbitrary symmetric matrix. In the case where 
all control variables have a similar ability to affect the 
trajectory, W is taken equal to the unit matrix, and DP^ be- 
comes the integrated square of the control variable perturbations 
6a (t). It might be noted that if Equation 5 is to have meaning, 
it is essential that all control variables have the same dimen- 
sions. To meet this condition, the control variables can be 
expressed in non-dimensional form. 

The constraint on control variable perturbation size repre- 
sented by Equation (5) is an essential element of the steepest 
descent process; for the optimum perturbation will be found by 
local linearization of the non-linear trajectory equations about 
the nominal path. To insure validity of the linearized approx- 
imation, the analysis must be limited to small control variable 
perturbations by means of Equation (5) which provides an integ- 
ral measure of the local perturbation magnitudes. 


Single Stage Analysis 

The steepest descent process has been outlined above. To 
implement this method, an analysis of all perturbations about 
the nominal trajectory must be undertaken. In the present 
report, all perturbations will be linearized; only first 
order perturbations in the control and state variables will be 
considered. The objective of the linearized analysis is 
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determination of the optimum control variable perturbation in 
the sense discussed in the previous section. 

Denoting variables on the nominal trajectory by a bar 


|a^(t) nominal » 
Xjj(t) nominal ■ 


where there are M control variables and N state variables. 

Now consider a small perturbation to the control variable 
history ,6a (t) ; this in turn will cause a small perturbation in 
the state variable history, 6x(t). The new values of the vari- 
ables become 

|a(t)j « jo(t)| -1^ jaa(t)! (8) 


x(t)f e ] x(t)[ ]ax(t) 


The nominal state variable and perturbed state variable 
histories can also be written as 

|x(t)| - |x(to)| + jf^x (t), a (t), t^l dt ^10- 

{x(t)| = J jf^x 6x, a -^Sa, t^j dt 

Subtracting Equation (10) from Equation (11) and using 
Taylor's expansion to first order, 

o * / I 


where 


f - f ^x(t), a (t), 


and where the repeated index indicates a summation over all 
possible values. Differentiation leads to 






or in matrix form 


[f]|4x| +[o]{«a} 


(14b) 


where 


, Mi ^ o Mi 

ij • 34 ^ 

th ^ T 4 4 M ^ M 4 th 


(15) 


Here the (i,j)^“ ^sjlement lies in the i^^ row and column of 
the matrices; F is an N x N matrix and G is an N x M matrix. 

The effect of these perturbations on pay-off, cut-off, and 
constraint functions must now be determined. A general method 
for obtaining these effects, known as the 'adjoint method,' 
Reference 13, is to define a new set of variables by the equations 


'x(t)]» -[K(t)j'[x(t)j 


(16) 


By specifying various boundary conditions on the X, the 
changes in all functions of interest can be found in turn. To 
show this pre-multiply Equation (14) by X' and Equation (16) 
by 6x' , transpose the second of these equations and sum with 
the first giving 


[xj.l4^(ax)|.[i]'j4x| 


p1 


6x] 


which may be written as 


, jit •['‘]'[°]W 


(17) 


(18) 


Integrating Equation (18) over the trajectory 


{x'4x}t - (x'«x}t„ ' jj'‘] '[“]{*“} 


dt 


(19) 


Now define three distinct sets of X functions by applying 
the following boundary conditions at t = T; 


d<t> 

axi 


{x(T)| 

[>'«]■ K 


- 

» 

= {Xq(t)} 

. [.,(«] 


(20a) 

(20b) 

(20c) 
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Equation (16) may now be integrated in the reverse 
(i.e., from T to t_) to obtain the functions, {XA(t)), 
and {X,^ (t) }. 

direction 
{ Xq ( t) } , 

Siibstituting each of these functions into Equation 
turn and noting that 

(19) in 


(21a) 


(21b) 


(21c) 

It follows that 



(22 a) 

«OfT -f 

*^o 

(72h) 

+|x^(tQ)J |ix(tQ)| 



(22c) 


Now, Equations (22) give the changes in pay-off function, 
cut-off function and constraint functions at the terminal time 
of the nominal trajectory; however. On the perturbed trajectory, 
the cut-off will usually occur at some perturbed time, T +7.T. 

In this case, the total change in the above quantities becomes 

aa - y + liCr) at 

I • 

v][.h dt + [x*(to)] 

11 



(23b) 

(23c) 


Equations (23) supply the change in pay-off, cut-off, and 
constraint functions on the perturbed trajectory. 


The time perturbation in Equations (23a) and (23c) may be 
eliminated by noting that, by definition of the cut-off function. 
Equation (23b) must be zero. 



(24) 


Substituting Equation C4) into Equations (23a) and (23c^ 


where 




(25a) 


(25b) 


(26a) 

(26b) 


Equations (25) reveal the significance of the X functions, 
originally defined by Equations (16) and (20). At time t^^, 

gives the sensitivity of <j>(T) to small perturbations in the 
state variables at to* Similarly, X({)f 2 (t) measures the sensi- 
tivity of (T) to small perturbations in the state variables at 
any time t. The sensitivity of the constraints dii^ to small 
state variable perturbations at any time is likewise defined 
by each row of the function 

A measure of the sensitivity of a trajectory to control 
variable perturbations can be obtained from the quantities X^'j^G 
and Xijj'j^G. Consider a pulse control variable perturbation at 
time t', that is, 6 (t-f ) , where 6 is the Dirac delta function. 
With this type of control variable perturbation/ it can be seen 
from Equations (25) that the changes in pay-off and constraint 
functions will be X^q (t ' ) ' G (t ' ) and X,j;j^ (t ' ) ' G (t ' ) , respectively, 
for fixed initial conditions. 

In order to apply the steepest-descent process, the perfor- 
mance function change, Equation (20a), must be maximized; subject 
to specified changes in the constraints. Equation (25b); and a 
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given size perturbation to the control variables, Equation (5). 
This can be achieved by constructing an augmented function in 
the manner of Lagrange which is to be maximized instead of dtji. 
For the present problem, the augmented function is 


” K H * [w*o)J { ix(to)^ 
*'^0 


(27) 


where the v are P undetermined Lagrangian multipliers, and y 
is a single undetermined Lagrangian multiplier. The objective 
now is to find that variation of the control variable history 
which will maximize U. 


Consider a variation of 6a , that is a 6 (6a). Then, 
it is always possible to write any 6a distribution in the form 

|5a|« I A(t)| k, or 

where A(t) prescribes the perturbation shape; and k, its mag- 
nitude. Now that part of Equation (27) which depends on 6a, 
the perturbation in the control variable, can be written in 
the form 

So that 


au 

97 


- j* L»*nj[o]{A(t)) dt [o]{A(t)} dt 


-T r 

+ 2ki [a( 1)J [wl |A(t)| dt 

IS 


( 30 ) 


or 


5U 


'to 

T 


+ 2M[_k . A(t)J [w]|ak . A(t)|| dt 

L^-ahJH ^ JM'[«]* ^'‘L‘“J[“] f >j ** 


(31) 


where it has been noted from Equation (28) that 

5 (6a) === A(t) 6k (32) 

Now, since Equation (31) holds for any A(t) 2 . it follows 
that it is a general relationship. Further, for U to be an ex- 
tremal, 6U must be zero. 


If 0 has been_maximized by means of a control variable 
perturbation 6a, 6U must be stationary for all small pertur- 
bations to the 6a, that is, for all 6 (6a) . The only way in 
which Equation (31) can be zero for all 6 (6a) is for the 
coefficient of 6 (6a) to be identically zero. That this last 
statement is true follows from considering the case where, over 
some finite time interval between t^ and T, the coefficient of 
6 (6a) is, say, positive. If this were the case, we could 
choose a 6 (6a) distribution that was also positive in this 
same interval and zero elsewhere between tQ and T. It would 
follow that U was also positive, and, hence, U could not be 
maximum. A similar argument holds when 6 (6a) is negative over 
any interval in to to T.; Hence, the coefficient of 6 (6a) must 
be identically zero in the whole interval to <. t <. T. This 
argument is essentially based on that presented by Goldstein, 
Reference 14. it follows that 


[M * l"J[ J[“]- 


(33) 
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Transposing, noting that W is synunetric, and solving for 6a, 


Substituting Equation (34) into Equation C25b) 
where 

■t 

and 

**0 

For subsequent use define the integral 

o 

The multipliers v can be expressed in terms of the 
y by Equation (35a) 

Substituting Equation (34) into Equation (5) 

^ ■ ii? { ” } n " J i * L - J 

Transposing the second term in the right hand side 

H * I " J 

Substituting Equation (37) in Equation (39) 


(34) 

(35a) 

(35b) 

(36a) 

(36b) 

(36c) 

multipliers 

(37) 

(38) 

brac)cet 

(39) 
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and noting that symmetrical gives 

um® - [i**j { v} * L"'*J {«} 


So that 


o.. - [^4'4T 

"" -Vopa - L4^j [1^,3-i i«| 


( 41 ) 


Substituting Equation (41) into Equation (37), the remaining 
Lagrangian multipliers are obtained in the form 


■ I -1 ■ -M‘ |W4*J iy : Mj 


(42) 


The optimum control perturbation is found by substituting 
Equations (41) and (42) back into Equation (34) and is 

f ^ 

+ [w]-i[o][x^n][i*J‘^{««} (43) 


With this equation the steepest-descent control pertur- 
bation has been determined. Perturbing the control variables 
according to Equation (43) gives the optimum change in the 
trajectory as discussed in the section entitled, "Problem 
Statement," with the added effect of changes in the initial 
value of the state variables included through the term in dg. 
The appropriate sign to use on the first term of equation (43) 
can be determined by evaluating d4>. Substituting the optimum 
control perturbation into Equation (25a) results in the 
equation shown on the following page. 
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- L^**J [ ^4''- { W) - [^o\ {«} 

+ [l^*J[l^^]-’-{ds} +K$)(»o>J { (44) 

As the quantity in the radical must be positive to assure 
the change in <p is real, it follows that the negative sign 
must be taken when minimizing the payoff function and the 
positive sign when maximizing the payoff function. 


An Alternative Analysis Using the Independent Variable 

for Cut-Off 

In the analysis of the previous section, it is implied 
that ahy function of the form 

n(xa(T), T) - 0 (45) 

will suffice to terminate the trajectory. While this is true 
in an analytic sense, in practice any function passing through 
zero more than once in the cut-off region may be difficult to 
employ for cut-off purposes. 


Nominal 


Perturbed 




Figure 1. — Double Valued Cut-Off Function 


Figure 1 presents a nominal cut-off function history which 
decreases monotonically . The perturbed cut-off function 
history, shown dotted, behaves in a different manner in that 
it passes through zero twice in the cut-off region. As the 
trajectory must be integrated numerically, there is a danger 
that cut-off will occur the first time passes through zero 
instead of the second, thereby introducing both errors in 
the linearized perturbations and preventing the build-up of 
the anticipated cut-off function history. 
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- LH [v]'" (*'*}) 

♦ [I**J[I,^J■^ {<»<’} ♦Kn<‘o)J { (44) 

As the quantity in the radical must be positive to assure 
the change in ^ is real, it follows that the negative sign 
must be taken when minimizing the payoff function and the 
positive sign when maximizing the payoff function. 


POINT MASS TRAJECTORY EQUATIONS 


Basic State Variables 

Preceding portions of this report derived a successive 
approximation scheme for computation of optimum trajectories 
generated by a set of first order differential equations. The 
analysis is quite general and holds for trajectories generated 
by any set of first order differential equations. The object 
of the present section is to specialize the analysis to point 
mass vehicle trajectory problems. This will be accomplished 
when a suitable set of state variables, together with their 
derivatives# the control variables, and the forces associated 
with the control variables has been specified. First, a suit- 
able coordinate system is selected, and Newton's Laws in this 
system are used to define the vehicle's motion. 

Several suitable coordinate systems are available for 
point mass trajectory computations. The basic set of coordi- 
nates used in the present analysis is a rectangular set rota- 
ting with the earth, (X^, Yq, Zg) . This coordinate system is 
illustrated in Figure 4. 


I 



Figure 4. — Basic Coordinate System 




The Xe and Ye axes lie in the equatorial piano, ih«j po.oit;ive 
Xq axis being initially chosen as the inters', ction of Uii; plane 
with the vehicle longitudinal plane at t = t^ . Yo i to 
the west of Xe# and Zq is positive through the South Polo. 
Denoting the radius vector from the center of the earth to the 
point mass vehicle by R / its magnitude is given by. 


1 * 1 - 

(45) 


The angle between R and the North Pole is given by 

<t>' 90 - <t>L (46) 

where tpjj is the latitude of the vehicle. As a result of the 
earth's rotation, an observer in the (Xe» Ye» Ze) system would 
detect an apparent motion of the point mass if it were at rest 
in inertial space. In time At the apparent displacement of 
such a vehicle would be 


‘’'apparent = « si" • “p 'it 
to the west. In vector notation 


(47) 


^^ajpparent “fcxwpAt = -<jj) x Rdt (48) 

This apparent displacement is independent of the v hide's 
motion and exists whether or not the vehicle is at rest in 
inertial space. In general, then, to an observer in the rota- 
ting coordinate system. 


(iR)e “ (*K)ln*rtlal + ( ^ R)app«r«nt (49) 

.*• Inertial " (^®)e (50) 

Dividing Equation (50) by AT and taking the limit, it 
follows that 


^Inertial “ (dt)e 
Ylnertial - 


(51) 

(52) 
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The vector • in Equation (51) / could equally well be taken 
as any vector? the arguments of Equations (45), ' to (52j ‘ still 
hold. Therefore, in general, for any vector quantity the oper- 
ational’ equality 




(55) 


can be defined. Applying Equation (142) to Equation (141), the 
inertial acceleration is given by 



(54) 


Now Newton's Law applies in inertial space so that in the 
rotating system 


I 

ffl 






(55) 


Here F is the total force acting on the vehicle. This vector 
equation can be expressed in component form using the relation- 
ships 


R “ X^.I + + Vk 


(56) 


«p ■ ■“'p’k 

F - F .1 + F .j + F, .k 
e ' e 


(57) 

(58) 


Here k are unit vectors aligned along the X^, Y^, and 

Zq axes, respectively. Equating components on either side of 
Equation (55) 
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- Xe * a-pY, - Xe 


- Ye - a^Xe - <-p^Ye 


e « 


These equations are not yet in a suitable form for the 
steepest descent analysis to be applied, for they are not in 
first order form. The transformation of Equations (59) to 
.(61) ' into first order form is immediately accomplished if 
the following quantities are defined as state variables : 


{*} 'lui 


where 


e u-.l + v«.i + w*.k 

ar- * * 


With this set of state variables the following expressions for 
the state variable derivatives are obtained from Equations (59) 
to (63.) . 

Xe - ue (64) 

■ 

*= (65) 

( 66 ), 

^ . 2o,pVe 4 0^^ X, (67) 
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These equations are In the same form as Equation (1) 
provided the total force is a function of the state variables, 
a set of control variables, and time. When the mass is vari- 
able, it too must be introduced as a state variable. Any 
expression for the rate of change of mass of the form 

■ - * ^ 

may be used in the analysis. The above state variables, Xe, 
Yg, Zq, Uq, Vg, We, and m will be referred to as the basic 
state variables. In certain problems it becomes necessary to 
specify additional state variables; these will be discussed 
later in this section of the report. 


Control Variables 

The total force acting on the vehicle has three distinct 
sources: first, aerodynamic force as a result of interaction 

between the vehicle surfaces and the planetary atmosphere; 
second, gravitational force as a result of vehicle and planetary 
mass interaction; and finally, thrust forces introduced by the 
vehicle propulsion system. 

Before aerodynamic forces can be computed , the atmospheric 
properties, vehicle velocity relative to the atmosphere, and 
vehicle attitude must be specified. Atmospheric properties can 
usually be specified as a function of altitude which, in turn, 
is a function of the state variables Xg, Ye, Z». Vehicle velo- 
city relative to the atmosphere is also a function of the state 
variables, for Ue, Ve, and w^ are the vehicle velocity compo- 
nents in a rotating system. The first and second factors 
determining aerodynamic forces are, therefore, functions of 
the basic state variables. 

The remaining factor entering into aerodynamic force 
determination, the vehicle attitude, is clearly not a function 
of the basic state variables. For, given the vehicle's position 
cind velocity, we are still quite free to specify its angular 
orientation in space. The angles which determine vehicle 
orientation may, therefore, be utilized as control variables by 
which aerodynamic forces may be modulated. Any set of three 
independent angles could be utilized for this purpose. Conven- 
tion suggests use of the vehicle angle-of-attack and angle-of- 
sideslip to orient the vehicle reference axis with respect to 
the velocity vector. Angle-of-attack, (a) , is the angle 

between the velocity vector and the vehicle reference axis when 
viewed in the vehicle side elevation. That is in a rectangular 
coordinate system, x, y, z with x along the vehicle reference 
axis, positive forward, y perpendicular to the vehicle plane of 
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symmetry, positive to starboard, and z completing a right hand 
system, a view normal to the x-z plane is considered. If u, v, 
w are the components of the vehicle velocity with respect to the 
atmosphere in this body axis system 

a - tan"^(g) '*71) 



Sideslip fmgle (0) is the angle between the velocity 
vector and the reference axis when looking down on the vehicle 
planform, that is along the z axis. In this case, 

0 - tan“^ (^) (72) 



Angle of attack and sideslip completely define the attitude 
of the vehicle with respect to the velocity vector# The third 
angle required to establish vehicle orientation in space is a 
rotation about the velocity vector. This last angle, bank angle 
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(B^) , will be. taken as zero when the vehicle plane of symmetry 
is vertical and the vehicle upright. Positive bank angle will 
be taken as a positive rotation about the velocity vector, as 
in Figure 7. 


Vertical plane containing velocity 


Vehicle plane of symmetry 




is the angle between the . 
aircraft z axis and the ^ 
vertical when viewed along 
the velocity vector 


Figure 7.- 


Bank Angle 


With the above set of angles to describe vehicle attitude, 
the velocity vector known and a given atmosphere, the aerodynamic 
forces can be completely specified. 

Returning to the second source of vehicle force, gravitation, 
from Newton ' s Laws , this is merely a function of position and 
mass. It is, therefore, completely defined in terms of the state 
variables and, hence, introduces no new control variable. 

The final source of vehicle force, thrust from the propul- 
sion system, involves the atmospheric properties, either due to 
the atmospheric back pressure degrading the vacuum thrust or by 
virtue of the atmospheric fluid used in the combustion process 
Which creates thrust. The propulsion unit efficiency may be 
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affected by Mach number and, hence, velocity so that thrust 
forces depend on the basic state variables of position and 
velocity in a similar manner to aerodynamic forces. If the 
propulsion system has a fixed orientation within the vehicle, 
the control variables introduced to describe aerodynamic forces 
suffice to describe thrust forces also. It may be, however, 
that the propulsion unit has a variable orientation within the 
vehicle. In this case, additional control variables to describe 
the relative position of the propulsion unit with respect to 
the vehicle are required. With vehicle attitude already spec- 
ified by a, 3 and two additional angles are sufficient 

to orient the thrust. These may conveniently be taken as the 
cone angle from the reference axis, inclination 

about the reference axis, 4>ip. This latter angle will be 
measured positively about the reference axis and be zero when 
the thrust force is perpendicular to the port side of the veh- 
icle plane of symmetry, as Illustrated in Figure 8. 




Figure 8. — Thrust Angles 

One other control variable for thrust remains to be speci- 
fied; this is the throttle setting, N, which serves to determine 
the propulsion unit power setting on variable thrust engines. 

In all then, to specify the forces acting on a point mass 
vehicle with a single propulsion unit, six control variables, 

Of Ba» <I>T/ Sind N are required. If there is more than 
one independently controllable propulsion unit, additional X<p, 
and N must be defined. 
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Coordinates and Coordinate Transformations 

Local geocentric-horizon coordinates. — Components of 
the planet-refetenced acceleration are irTtearated to obtain 
the planet-referenced velocity components ( xq , 2e) • Vehicle 

position in this coordinate system is determined by integration 
of these velocities. Vehicle position in the planet-referenced 
spherical coordinate system will now be determined. The spher- 
ical coordinates are longitude, geocentric latitude, and dis- 
tance from the center of the planet. Angle "C" represents the 
change in vehicle longitude and may be written 

C « 0T - 6t 
Lq l 

Angle C is related to the vehicle position by the expression 

-1 Ye 

C = Tan (rr-) ^ r74l 


The relationships are illustrated in Figure 9. 


Local-geocienlric 
horizon coordinates 



Figure 9. — Relation between Local-Geocentric, Inertial and 
Earth-Referenced Coordinates for Point-Mass Problems 


To describe body motion relative to the planet, a local- 
geocentric-horizon coordinate system is employed. The Zg-axis 
of this system is along a radial line passing through the body 
center of gravity and is positive toward the center planet. 

The Xg-axis of this system is normal to the Zg-axis and is posi* 
tive northward; Yg forms a right-handed system. Figure 9 shows 
the relation of this coordinate system to the other systems 
employed. 



To locate the Xg-Yg-Zg axes with respect to the Xe-Ye-Ze axes, 
first rotate about Zq by an angle (180® + C) and then rotate 
about Yg through the angle (90® - (()l) • The first rotation 
defines the intermediate coordinate system shown in Figure 10. 
The transformation is given by 


ix' 


Coa (100° ♦ C) Sin (l80° + C) 0 




B 

-Sin (100° * C) Cc5 (100° + C) 0 





0 0 1 




or 


lx’ 



m 




-Cos 

Sin 

0 


C 

C 


-Sin C 0 

-Cos C 0 
0 1 



(75) 


The second rotation is shown in Figure H. The transfor- 
mation matrix for the second rotation is given by 


or 


lx 


Cob ( 90 ° - ^l) 0 -Sin ( 90 ° 7 


m 

.0 1 ■ 

0 




Sin (90° - 0 Cos (90° - 

Iv 


Sin 4 l C "Cos 


ix' 

> 

u 

0 10 


> 

K 


Coa C Sin 


^Ze 


ix' 


(76) 



Figure lo- Intermediate Coordinate System Transformation 
from Earth Referenced to Lochi Geocentric Coordinates 





T 


Z 

e 

Figure ii- Final Rotation in Transformation from 
Eartti-Referenced to Local-Geocentric Coordinates 
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In this analysis, a positive rotation is defined in the 
sense adopted for vector cross products in a right-handed system. 
That is^ a positive rotation about the z-axis occurs when the 
x-axis rotates into the y-axis; positive rotation about the 
x-axis when y-axis rotates into the z-axis; and positive rota- 
tion about the y-axis when the z-axis rotates into the x-axis. 

The intermediate coordinate system (X', Yg, z^) is eliminated 
by the method of successive rotation. The complete transfor- 
mation is given by 


1y 


Sin ♦l 0 -Cos 


-Cos C -Sin C 0 


1y 

■''■e 

> 

= 

0 1 0 


Sin C -Cos C 0 





Cos 0 Sin 4*^ 


0 0 1 




C7I) 


This can be reduced to the single transformation matrix 


ix 


^g 

II 

iz 

^g 

1 


-Sin ^ -Sin Sin C -Cos 

Sin C -Cos C 0 

-Cos Cos C -Cos Sin C Sin 4>jj 


lx. 




(78) 


which defines a direction cosine set (i, j, k) by the equation 


1y 

^g 


ii h 


IXe 

ly 

ig 


i2 J2 ^2 



I7 


is ^3 




(79) 


Planet referenced velocity in the local-geocentric coordinate 
system is given by 


and 


ig 


ii 

h 





Yg 

= 

i2 

J 2 

k2 


Ye 


Zg 


h 

•J 3 

k3 


Ze 

(80) 

^g 

m 

’ g 

g 


• 

(81) 


Flight path angles are computed by 

0 ■ tan" 




(82) 
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and 


y ■ 


a in 


-1 



( 83 ) 


Here a is the heading angle and \ is the flight path angle. 

Wind Axis Coordinates . — Aerodynamic and thrust forces 
for point-mass problems are conveniently summed in a wind-axis 
coordinate system, Ya, Za) . Since the equations of motion 

are solved in {Xq, Yg, Zg) coordinates, the wind-axis components 
of force must then be resolved into this basic system. 

When winds exist, defined by atmospheric velocity compo- 
nents along the local geocentric axes, vehicle velocity rela- 
tive to the atmosphere is the vector difference of vehicle 
geocentric velocity and wind velocity. The wind axis system 
is then determined by the vehicle airspeed, and the flight 

path angles relative to the atmosphere Xa wind 

velocity is zero, Va = Vg, Xa = X and oa = o. If there is 
a wind, with velocity components (Xg^,, Yg^, Zg^) , then 


Va-\ 

1 <VV>' * * ' VV'" 

(84) 

Ya = 

sln-1 |-(Xg-X^p/vJ 

(85) 

'"a ' 

tan"^ 1 (Y -Y )/(X -X ) 

L g gw ' g gw'J 

(86) 


Forces are first resolved from v/ind axes to the local 
geocentric coordinates. The wind axes are defined relative 
to the local geocentric axes by three angles: heading, qa? 

flight path attitude, Xa/ defined above; together with angle, 
Ba. 


h 


Figure 12. — 

Relationship between Local-Geocentric Axes and Wind Axes 
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The resolution of forces from wind axes to local geocentric then 
becomes 
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( 91 ) 


For the rotating-planet, the local geocentric components must 
be resolved into the Xg-Yg-Ze system. The required direction 
cosines are given by Equation 
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( 92 ) 


The combined transformation from wind axes to local geocentric 
can be defined as a single matrix transformation 
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Body-axis coordinates. — Origin of this system is the ve- 
hicle center of gravity with x-axis along the geometric longi- 
tudinal axis of the body. Positive direction of the x-axis is 
from center of gravity to the front of the body. The y-axis is 
positive to starboard extending from the center of gravity in 
a water-line plane. The z-axis forms a right-handed orthogonal 
system. To permit the use of body (x, y, z) axes aerodynamic 
data, and to convert the body axes components of thrust to the 
wind axes system, a coordinate transformation must be made. The 
coordinate transformation shown in Figure 13 involves rotation 
first through angle of attack, a, then through an auxiliary 
angle, 0'. Noting that 


tan 0 ' - 


V 

u 


cos a ■ tan 0 cos a 


(94) 
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the transformation is 


X* cos a 0 sin a x 

y» ■ 0 1 0 y 

z' -sin a 0 cos a z 

cos 0' sin 0* 0 x’ 

** -sin 0’ cos 0' 0 y' 

Za 0 0 1 z ’ 

■ cos. 0' cos a sin 0' cos 0' sin a x 

-sin 0' COB q cos 0' -sin 0' sin o y 

-sin a 0 cos a z 

which defines the (u, v, w) direction cosines 
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The relationship between body and wind-axes aerodynamic coef- 
ficients is now established. Note the negative directions of 
the coefficients relative to the axes. 

Inertial coordinates . — The selected inertial coordinates 
coincide with the earth references system at time 

zero. At a later time they differ by the rotation of the earth 


34 


ojpt. The transformation between inertial velocities and planet 
referenced velocities is derived as follows. 

Let R be the displacement of the point-mass, {See Figure 9). 

In inertial coordinates 

R « Xlx + Yly + ZI2 (98) 


and 


V - R = ilx + Yly + ziz c^9) 

In planet-referenced coordinates 

« ■ 5teix„ + + 2elz 

^ e e "e 

However, due to the rotation of the Xg, Yg, Zg coordinate system, 
the velocity is 

— £R — — 

V-R-^+OJpXR (100) 

where 

* ^e^Xg + ^e^Yg + Ze^Zg (101) 


The planet's rotation is about the Z-axis which is also the Zg- 
axis. Therefore, 


Wp ® “t^plz = “‘‘‘p^Zg 
and the required cross product is 





= (YeWp)ixe ~ 


( 102 ) 


If Equations , (99), (101) , and (102) are substituted into 

Equation (189) , it / follows that 


Xix + Yly + ZI 2 « (Xg + a3pYg)Ixg + (Yg - a)pXe)ly^ + (Ze)lZg 

(103) 

The relation between the unit vectors in the inertial system 
and unit vectors in the planet referenced system are obtained 
by a single rotation about the Z-axis. 
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The transformation matrix is 



( 104 ) 


The transformation from planet-referenced velocities to inertial 
velocities is made with the inverse of the matrix of Equation 
C104) f and the component relations derived in Equation ■(103) 


• 
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Cos (dpt Sin (dpt 0 


Xg ♦ 
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-Sin (dpt Cos (dpt 0 


Ye - “pXe 
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0 0 1 


Ze 


(105) 


The components of inertial velocities are used to calculate the 
inertial speed of the body as 


V. 


2 *2 '2 
X + Y + Z 


C106) 


Equation (106)^. is valid regardless of the inertial coordinate 
system involved. 

Local-geocentric to geodetic coordinates . — Positions on 
the planet are specified in terms of geodetic latitude and 
altitude (for a given longitude) while the motion of the body 
is computed in a planetocentric system which is independent of 
the surface. In the computer program, flight-path angle X 
and heading angle a are calculated with respect to the local 
geocentric coordinates. By definition Xp and oq are angles 
measured with respect to the local geodetic. Although the 
maximum difference that can exist between the two coordinate 
systems is 11 minutes of arc, it may be desirable to know X^, 
and OD more accurately than is obtained when measured from 
the local geocentric. 

It will be necessary to resolve the geocentric latitude 
to geodetic latitude for an accurate determination of position. 
Figure 14 presents the geometry required for describing the 
position of a point in a meridian plane of a planet shaped in 
the form of an oblate spheroid* 
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Figure 14, — Planet-Oblateness Effect on Latitude and Altitude 

It is apparent from this figure that the most significant 
difference between the geocentric referenced position and the 
geodetic position is the distance AB on the surface of the 
reference spheroid. The distance can be defined by a knowledge 
of the angle the geocentric latitude; (j>g, the geodetic 

latitude; the corresponding radii; and the distance OC. 

The relationship between the geocentric and geodetic lati- 
tude of a point on the surface of a planet which is an oblate 

The equation for the surface 


= 1 

(107) 

The tangent of the geodetic latitude can be found by determining 
the negative reciprocal of the slope of a tangent to this ellipse 
The expression for this tangent is 

- 31 ^) 

dX 

Note that is a negative number in the northern hemisphere. 


Re Zb 


V 


(108) 




spheroid is obtained as follows, 
in a meridian plane is 




i- 


. I 


The tangent of the geocentric latitude of point B 

T*" ♦l, ■ - 
. ® . 3<B 


IS 
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Substituting Equation 
required relation 


Tan ♦g 


(109) into Equation (108), gives the 



Tan ♦Lg 


■(110) 


The expression for the radius of the planet at point B in 
terms of the geocentric latitude of the point and the equato- 
rial and polar radii is obtained by the rectangular to polar 
coordinate transformation 



(111) 

Xb = R^Lg ^Lg 

, (112) 


and, solving for Ra by substituting Equations (111)/ and (112)^ 
into Equation (^io75»^’^9 gives 


“ V Rp2 Sin^ + Re^ Cos^ 4^ 
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(313) 


cos 

cos 




i./vi 
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It may be seen from Figure 14 that 




ITT’'' * OF' - OB' 


(114) 

or 

h sin ^g * OP sin “ ^♦Lg 

♦tg 

^i:ii5) 

Likewise 

B'rp” s CF' - 00" 


(116) 

or 

h cos ♦g = OF cos cos 


(117) 

If Equation ■ 

(115) is divided by Equation 

(117)/ 

and then the 


quotient is divided by tan /there results 
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o g 6 ’» 

"V *L,)/ 

(^) (tan ♦j^/tan ♦,) - (|) * [l - (Ep/s,)^] (^^^) («* /op) j^O) 

Let 

U ■ (Rg tan /Rp tan 

* (Rp tan i|>g/Re tan (j»j^) (121) 

Then it follows from Equations (113) ' and (120)^? that 

U » (jj^) + [E^/OT] [U/ ^ sin® * cos® H - (Rp/R^> 1 

Equation (122)- is solved by an iterative scheme. 

Then 


, R U 

* tan"-^ [(-| — ) tan 4.^^] 
P 


(123) 


The flight-path and heading angles corrected to the local 
geodetic latitude are computed by 


.-1 =sin-l (:...^6 -|y*fi-*L)}j 


(124) 


Since the magnitude of vector is equal to the magnitude of 
vector V, ^ 


9l 


and 


Ojj = Sin 




Zg(*g - 


( 125 ) 
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Auxiliary Computations 


In addition to the computations which can be made from the 
problem formulation as presented in preceding sections, several 
other quantities are available as optional calculations. 


a. 

b. 

c . 

d. 

e . 

f . 


Planet-surface referenced range, R, 
Great-circle range, R 


D 


Down- and cross-range, Xq and Yq 
T heoretical burnout velocity, V^heo 
Velocity losses, Vp, Vp, and Vml 

Orbital Variables and satellite target 


Planet-surfaced referenced range . ?»- The total distance 
traveled over the surface of the planet is computed as the 
integrated surface range. If the distance traveled by the 
vehicle over a given portion of the trajectory is 




/ 


,t2 


Vg dt 
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then the curvilinear planet surface referenced range is 

■ .tg 


■/ 




L Vg Cos Y dt 


(127) 


The flight-path angle, X, is referenced to local geocentric 
coordinates for this computation. 

Great-circle range . — Great-circle distance from the launch 
point to the instantaneous vehicle position, Rg, may also be 
required. Expressions for this distance are derived as follows. 

By spherical trigonometry, (see Figure 15 j 
Cob « Cos (90 -^l)Cos(90-^l^) + Sin(90-*L)Sin(90-«>Lo) Cos(eL-eL^) ' 


or simplifying 
R 


Cos ■ Sin Sin ♦lq 


(128) 

(129) 


Therefore, 

Rg ■ R* Cos“^ j^Sin Sin + Cos Cos 4^0 Cos (6l“®Lo)] ‘ (130) 
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Launch point 




ALOl. 


Equator 


Figure 15.- — Great-Circle Range 


However, since the planets are generally oblate spheroids, R* 
is not a constant radius. An approximation may be obtained by 
averaging the planet's radius at the launch point and at the 
vehicle's position. Therefore, define the average radius, R' , 




(131) 


and the surface-referenced great-circle range from the launch 
point to the vehicle is 

Rg = ^ Cos"^ rSin4>L Sin Cos Cos Cos(0L-eL^)j 

Down- and cross-range . — Down- and cross-range from the 
initial great circle can be determined. The initial great 
circle is determined from the input quantities Oq, 4>Lo» and 



0L (see Figure 16) Then the cross range of a particular 
trajectory point is defined as the perpendicular distance from 
the point to the initial great circle. The downrange is then 
the distance along the initial great circle from the initial 
point to the point P at which the cross range is measured. From 
the spherical triangle, Figure 16 , the great circle range LF to 
the point F is computed by Equation 


The right spherical triangle LPF is solved for the down- 
range, Xp, and the cross range, Yp. 


Xp ■ R’ cos 


-1 


(: 


Cos LF 


Cos (sin"^(3in LF sin 0 ) 
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Yn * R' sin”^ (sin LF sin C) 


(134) 


where 


5 « C - 0, 


R* is defined by Equation (131) 


(135) 



42 


i 


Theoretical burnout velocity and losses . — For trajectory 
and performance optimization studies/ it is convenient to know 
the theoretical burnout velocity possible and the velocity 
losses due to gravity, aerodynamic drag, and atmospheric back 
pressure upon the engine nozzle. These quantities may be com- 
puted as follows: 

Theoretical Velocity 


t2 


theo = 
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^VAC 

m 


dt 
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Speed Loss Due to Gravity 


grav 


*1 


g2g Sin Y dt 
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Speed Loss Due to Aerodynamic 
Drag 


Vd = 


/ 

^1 


D 

m 


dt 


Speed Loss Due to Atmosphere 
Back Pressure Upon the Engine Nozzle 


/ 


PA 


m 


dt 


(138) 


(139) 


'ML 


Maneuvering Losses 
'^VAC ■ 




(140) 


The resultant velocity Vg(t2> is obtained by adding the components 
computed to the initial value Vg(tj^) 


Vg(t2) - Vg(ti) + V,heo + Vav + Vq + Vp 


(141)' 


The maneuvering losses are valid only if X t is zero for the 
engine. 


43 


Orbital variables and satellite target . — Certain functions 
of use in orbital trajectory calculations have been added to the 
point mass equations of motion used in the Steepest Descent Opti- 
mization Program. These functions permit the specification of 
terminal conditions in inertial space when this is convenient. A 
further set of functions will permit rendezvous cklculations with 
a satellite in a circular orbit about a central planet. 


Orbital variable calculations commence immediately after 
the calculation of vehicle inertial velocity. Flight path 
angles in inertial space are computed from the expressions 



(142) 


(143) 


The inclination angle, i, is the angle between the plane con- 
taining the velocity vector and the center of the earth, and 
the equatorial plane. 


-Z 



Figure 17 .— Orbital Plane Geometry 
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Applying spherical trigonometry to Figure 17 , we obtain the re- 
lationship 


COB i ■ cos Bln Oj (144) 

The difference in longitude between the vehicle and the ascending 
node, V, is given by 

tan V = sin tan Oj 
The inertial longitude is given by 

®I ®L " 

and the inertial longitude of the ascending node by 

ft = 01 - V (147) 

It is convenient to know the central angle, u, in the orbital 
plane. Measuring from the ascending node, we obtain 

tan 6r 

tan u = ^ (148) 

cos Oj 


(145) 

(146) 


The orbital variable, calculation introduces positional and 
velocity information from a second body. This body is a satel- 
lite considered in a circular orbit about the earth. Its orbital 
height, hg, is specified and remains constant. Position in the 
orbit is computed from an initial central angle, ^is , by the 
expression ° 

= ♦sq (149) 


The satellite angular velocity is obtained from the satellite 

inertial velocity, , where 

s 



(150) 


where yg is the gravitational potential constant and Re the 
earth radius. It should be noted that Fquation (ISO) assumes 
a spherical earth; for the earth radius is taken as constant, 
and none of the higher order gravitational harmonics are in- 
cluded. Knowing Vcg, it follows that 


“s 




(151) 


The variables of this section provide sufficient information to 
either rendezvous with or terminate the trajectory in a speci- 
fied position relative to the satellite. 


VEHICLE CHARACTERISTICS 


Methods by which the aerodynamic, propulsive, and physical 
characteristics of a vehicle are introduced into the computer 
program are presented in this section. Form and preparation of 
the input data are discussed, together with methods by which 
stages and staging may be used to increase the effective data 
storage area allotted to a description of the vehicle's proper- 
ties. 


Aerodynamic Coefficients 

The primary objective of the aerodynamic data input sub- 
program is to provide for a complete accounting of the various 
contributions to the aerodynamic forces and moments regardless 
of the flight conditions of the vehicle being considered. Two 
techniques are available for use in the digital computer pro- 
gram: (a) an n-dimensional table look-up and interpolation and 

(b) an m-order polynomial function of n variables prepared by 
"curve fit" techniques. In the first method, the proper value 
for each term is obtained by an interpolation in "n" dimensions 
where the number of dimensions is taken to be the number of 
parameters to be varied independently plus the dependent vari- 
able. This method has the advantage of accurately describing 
most non-linear variations with reasonable preparation effort. 
The amount of storage space which must be allocated to such a 
method, however, can achieve unreasonable proportions and may 
require substantial computing time for the interpolation as the 
number of dimensions are increased. The second method has essen 
tially the opposite characteristics; that is, a large amount of 
data may be represented with a small amount of storage space, 
and computation time is held to reasonable limits, but the data 
variations which may be represented must be regular. A substan- 
tial amount of effort can be required for the preparation of 
data by a curve-fit technique. Both these methods are very con- 
venient when the amount of data to be handled is moderate, but 
tend to become unmanageable when large amounts of data are 
required. This usually occurs when the program, having several 
degrees of freedom, is committed to one or the other of these 
two techniques. Therefore, the computer program incorporates 
both of the techniques discussed as a compromise to take advan- 
tage of the more desirable features of both. To do this, a 
general set of data equations have been programmed which define 
each of the aerodynamic forces. In general, the coefficients 
for these equations will be obtained from a curve-read inter- 
polation. Several simplifications may be made to the equations 
depending on the flight condition and vehicle to be considered. 
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Often the particular application will not require some of 
the terms listed in order to describe the flight path and vehicle 
under consideration. The subprogram is arranged so that the com- 
puter will assign a constant value to any curve for which the 
data has not been supplied. For most curves, the constant value 
will be zero. This technique may be used to reduce the time 
required for the preparation of data. Values intermediate to 
those introduced in a tabular listing will be obtained by linear 
interpolation. 


Aerodynamic Forces 

Aerodynamic forces are customarily defined by three mutu- 
ally perpendicular forces. These are lift (L) , drag (D) , and 
side force (Y) . Lift force is perpendicular to the velocity 
vector in a vertical plane; drag force is measured along the 
velocity vector but in opposite direction; side force is meas- 
ured in the horizontal plane, positive toward the right, pro- 
vided the bank angle is zero. If the bank angle is not zero, 

L and Y will be rotated by -B;^ about the velocity vector. 
Coordinates are shown in Figure 18. 



Figure 18. — Aerodynamic Forces - Wind Axes 
These forces may be expressed in the form 


L = q(V,hJ SC£(V,h,a,8) 

C152) 

D = q(V,h) SCD(V,h,a,0) 

(153) 

Y = q(V,h) SCyCv.h.a.e) 

(154) 


where q is the dynamic pressure and S a convenient reference 
area. The aerodynamic coefficients C^, Cq, and Gy may be ex- 
pressed in terms of the aerodynamic derivatives. 
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a + Cl^ 2 “[a | + Cl^ | b| 


+ CLg2B2 + “|e| 

C155) 

|a|+ Cd^2 |b| 



(156) 

Cy = Cy^ + Cy^|“|+ Cy^2 0^ + 


+ Cyg2 B|B|+ Cy^ 3 |a|B 

(157) 


Alternatively, the aerodynamic derivatives may be expressed 
as tabular functions of Mach number (M^) , ot, and 3, that is, 
functions of the state variables and the control variables. 

On occasion, it may be convenient to measure the aerodynamic 
forces in the body axis coordinate system introduced in a pre- 
ceding section, pages 28* to 30.. in this case, normal force, (nf), 
is measured along the -z axis, side force (y) along the y axis, 
and axial force (a) along the -x axis, as in Figure 19. 



I t 


2 2 
Figure 19. — Aerodynamic Force in Body Axes 

The specification of forces in the body axis system is 
similar to that in the wind axis system 


ny 

» qSCN 

(158) 

a 

= qSCA 

(159) 

y 

■ qSCy 

(160) 
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where the body axis aerodynamic coefficients are 
Cn ■ CNq + ♦ Cn^2 a |a | 

^ e2 ^ a |e| 

Ca * CAq + Ca^|oi| + Ca^2 

+ CAg |b I + CAg2 6^ + 1“6| (162) 

Cy - Cy^ ♦ Cy__|ci| ♦ Cy^2a2 

* V ®M * “VaB 


Thrust and Fuel Plow Data 

The techniques employed to introduce thrust and fuel-flow 
data into the equations of motion are developed in an approach 
similar to that employed for aerodynamic data. An n-dimensional 
tabular listing and interpolation technique is used with the 
independent variables being defined by the type of propulsion 
unit being considered. For the present formulation y the propul- 
sion units are grouped into the following options: (1) rocket# 

(2) air breathing engine. 

Propulsion option (1) rocket . — The thrust of a rocket 
motor IS assumed variable with stage time, altitude , and , if the 
rocket is controllable, it will also vary with throttle setting. 
The altitude effect is determined by the exit area of the 
nozzle, Ag, and the ambient pressure, P. If the thrust is 
specified for some constant ambient air pressure, the altitude 
correction can be calculated within the subprogram. If the 
rocket motor is uncontrolled, the vacuum thrust, in pounds , will 
be introduced by a tabular listing as a function of time , in 
seconds, and corrected as follows: 

T = Max [T - PAg, 0] 
vac e 

The propellant consumption rate is specified by a tabular 
listing in slugs per second as a function of time, in seconds, 
for the single engine options, or computed from the thrust and 
the engine specific impulse, Isp/ for the multiple engine options 

If the rocket is controlled, the propellant mass flow rate nif 
is introduced by a tabular listing as a function of throttle 
setting. The thrust is then specified by a tabular listing as 
a function of mass flow rate. 
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Propulsion option (2) air breathing engines . — An air- 
fare atiHingen^Tne"^ strongly affected by the environmental 
conditions ur^der which it is operating. Engines which would 
be grouped in this classification are turbojets, reunjets, 
pulse jets, turboprops, and reciprocating machines. The param- 
eters considered significant in the program are 

(a) Altitude (h-ft) 

(b) Mach number (Mj^) 

(c) Angle of attack (a-degrees) , and 

(d) Throttle setting (N-units defined by problem) 

Both the thrust and fuel flow are functions of these variables. 
In order to accommodate these variables, a five-dimensional 
tabular listing and interpolation are used to obtain both thrust 
and fuel flow. The thrust has no further correction as the 
effects of all parameters are assumed included in the interpo- 
lated value. 

Engine perturbation factors . — The engine options include 
provision for two data scaling factors for use in parametric 
studies. These are in the form 


T = T^^^ + 


(165) 


Components of the thrust vector . — The equations used to 
reduce the thrust vector to its components along the body axes 
are 

T^ = T cosXip . (166) 


Ty = -T sinX^ cos<j>^ (167) 

and 

T^ = -T sinX^ sin(|>^ (168) 


and Xt are defined and explained in the control variable 
section. 


Reference weight and propellant consumed . — ^Rate of change 
of vehicle mass, m, is ^ set equal to the negative of the total 
mass flow rate, -m^. m is integrated to give variation of 
vehicle mass, m. The instantaneous mass is used in the compu- 
tation of the body motion. The reference weight is obtained 
by an auxiliary calculation 

- m(32.174) (169) 
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The propellant consumed is computed as 

mf ■ m^ - m (17-00* 

where is a reference mass input equal to the initial vehicle 
mass 


Stages and Staging 

A problem common in missile performance analyses and 
encountered frequently in airplane performance work is that 
of staging or the release of discrete masses from the contin- 
uing airframe. The effect of dropping a booster rocket or 
fuel tanks is often great enough to require that the complete 
set of aerodynamic data be changed. Configuration changes at 
constant weight, such as extending drag brakes or turning on 
afterburners, may also require revising the aerodynamic or 
physical characteristics of the vehicle. Another use of the 
staging technique is possible with the present computer pro- 
gram which does not involve physical changes to the configu- 
ration; this technique may be used tcj revise the aerodynamic 
descriptors as a function of aerodynamic attitude or Mach number. 
With this use of the stage concept, accurate descriptions of 
the forces acting upon the vehicle may be maintained over wide 
attitude ranges, if required. 
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VEHICLE ENVIRONMENT 


The models for simulating the environment in which a 
vehicle will operate are presented in this section. This 
environment includes the atmosphere properties, wind velocity, 
and the field associated with the planet over which the ve- 
hicle is moving. The shape of the planet and the conversion 
from geodetic to geocentric latitudes are also considered. 

In the discussions which follow, the descriptions of vehicle 
environment pertain to the planet Earth. The environmental 
simulation may be extended to any planet by replacing appro- 
priate constants in the describing equations. 


Atmosphere 

The concept of a model atmosphere was introduced many 
years ago, and over the years several models have been de- 
veloped. Reference 15/" outlines the historical background of 
the gradual evolution of the ARDC model. The original (1956) 
ARDC model (Reference 15*) was revised to reflect the density 
variation with altitude ' that was obtained from an analysis 
of artificial satellite orbit data. This revision is the 
widely used 1959 ARDC Model Atmosphere and is the basic option 
in the present program. 

The advantage of a model atmosphere is that it provides 
a common reference upon which performance calculations can 
be based. The model is not intended to be the "final word" 
on the properties of the atmosphere for a particular time 
and location. The atmosphere properties are quite variable 
and are affected by many parameters other than altitude. 

At the present time, the "state-of-the-art" is not advanced 
to the point where these parameters can be accounted for; 
it may be several years before the effects of some parameters 
can be evaluated. 

1959 ARDC Model Atmosphere .-- The 1959 ARDC Model Atmos- 
phere is specified in layers assuming either isothermal or 
linear temperature lapse-rate sections. This construction 
makes it very convenient to incorporate other atmospheres, 
either from specifications for design purposes or for other 
planets. The relations which mathematically specify the 
1959 ARDC Model Atmosphere are as follows (Reference 16 ) 
the 1959 ARDC Model Atmosphere is divided into 11 layers as 
noted in the table below. 
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Layer 

Hh -Lower Altitude 

Upper Altitude 


{Ceopotentlal) 

(Geopotentiai) 

• 

Meters 

Meters 

1 

. 0 

11,000 

2 

11,000 

25,000 

3 

25,000 

47,000 

4 

47,000 

53,000 

5 

53.000 

79,000 

6 

79,000 

' 90,000 

7 

90,000 

105,000 

3 

105,000 

160,000 

9 

160,000 

170,000 

10 

170,000 

200,000 

11 

200,000 

700,000 

For layers 1 , 

3 , 5 , 7 , 8, 9 , 10 , and 11 , 

a linear molecular 

scale temperature lapse-rate is assumed 
equations are used: 

and the following 


®g? 

m 

1 t 

.iU4oa 

• ,3045h/63567o6 

Meters 

(171) 


m 

(%>b * %(Hgp - Hb) 

m 

OR 

(172) 

T 

m 


A - B tan-1 ^ 


(173) 

P 

m 

Pb 

'l +Ki(Hgp - Hb)j -‘'2 


(174) 

P 

m 

■ 1 

[1 + Ki(Hgp - Hb) j-(l-«2) 

Slugs/Pt .3 

(175) 

Vs 

m 

49. 


Ft. /Sec. 

(176) 

¥ 

m 

2.269681 x 10-8 

Ft .2/Sec. 

(177) 


For the isothermal layers 2 , 4 , and 6, the following changes 
are made 

P . pj^e‘’'3<«eP-»b) (178) 
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Values of the temperature, pressure, density, and altitude at 
the base of each altitude layer are listed below along with 
the appropriate values K^, K 2 / and K^. 


kosc 


QiMntlty 

1 

2 3 


5 

6 


.. 223 $£ 9 X 3 X 10 ** 

* 0 .I38J165BO X 10-* 0 

-.15920187 X 10 * 

^ 0 


.$. 2 $ 6 U 22 

0 U. 38826} 

0 

-7,9921765 

0 

«3 

0 

, 1 }T 688}2 X 10*3 0 

.12086887 X 10*3 0 

,206231.1.2 X io ‘3 

IV. 

jiB.eaa 

?S‘I 968 389.968 

}06.t88 

506.788 

296.188 

Pb 

2 U 6.2170 

kTa. 67}99 } 1 . 97 }U 8 

2 . 51 }Ii 578 

1,2160383 

2.1082U85 X 10 "® 

“b 

2.37692 X 10 "^ 

7.0611076 X 10*'' 7.7663892 X 10'^ 2.8803201 X 10"^’ 1 . 39 >‘ 7 l 2 } X 10 '" 

k.ll 900 U 2 X le*"'’ 

■b 

0 

UOOO. 2 } 000 . 

67000. 

53000. 

79000. 



jar-v. 




7 

8 

9 

10 

,11 


. 2 lili. 5 aiil X 10 ' 

,86629910 X 10 *^ 

. 75 b 3 *H 23 X 10*3 

.350711.76 X 10*5 

.22212911. X 10*^ 


8.3bU966 

1.7062397 

3 .bl 61 i 79 b 

6.8329589 

9.7613696 


0 

0 

0 

0 

0 

»b 

296.188 

iio 6 .l 88 

2386.188 

2566.188 

2836.188 

Pb 

2 . 18 U 75 W X 10 ' 

■3 1.9564912 X 10 *^ 

7 ,} 60 b 667 X 10 *^ 

5,697l6Mi X lo"^ 

2.97697U6 X 10*^ 

Pb 

k.. 26 ll. 8 s 6 X 10 

•9 2 , 232 U. 2 k X 10 '^® 

1 . 8 b 588 li 9 X 10 **X 

1.3387990 X 10*^ 

6.1150607 X 10*13 

■b 

90000. 

10 } 000 . 

160000. 

170000. 

200000 . 


Values of the appropriate constants to be applied in the tem- 
perature equation. Equation , are listed below, 

Hgp(Ka) A B C D 

0-90 1. 0. 

9O-I8O .75951115 . 174 i 64 oU 220,000. 25,000. 

180-1200 .93578678 .27396592 180,000. 140,000. 

U, S. standard Atmosphere, 1962 . — The part of the U.S. 
Standard Atmosphere, 1962, below 90 kilometers geometric 
altitude (295,276 ft. altitude) is defined in the same way 
as the 1959 model — by the hydrostatic equation and a piecewise 
linear variation of temperature with geopotential altitude. 
Equations (171)' to ;(179>'‘ are , therefore, applicable with a 
different set of constants. These constants, based on the 
published tabulation of atmosphere properties (Reference 17 ) 
at the base altitudes, are presented below. The 1962 model 
uses a different set of relationships above 90 kilometers. 

These have not been included. The tables define 1962 model 
properties between sea level and 295,800 ft. geometric altitude. 
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Values of the temperature# pressure# density# and altitude 
at the base of each altitude layer are listed below along with 
the appropriate values of K2/ and K^* 


Layer 


Quantity 

1 

2 

3 

1* 

Kl 

-.2255877 X 10 “*^ 

0 

.1*80121*06 X 10’5 

.12199559 X 10 ' 

K 2 

-.5255871 X loi 

0 

.3281*1*801 X 1Q2 

.122021*70 X 10 ‘ 

K3 

0 

.1576958 X 10-3 

0 

0 

% 

518.67 

389.97 

389.97 

1*13.101* 

I'b 

2116.217 

1*72.6812 

111 *. 31*31 

17.22518 

Pb 

.2377002 X 10"2 

.7061512 X 10-3 

.1708202 X 10-3 

.21*29209 X 10 “' 

Hb 

0 

10999.1*71* 

19999.191 

3235l*.65l* 



Layer 


6 

Quantity 

5 

6 

7 


0 

-.7383899 X 10“5 

-.1572230 X 10 -^ 

0 

Kg 

0 

-.1709562 X lO'*’^ 

8602817 X 10 

0 

■ ^3 

.1262323 X 10“3 

0 

0 

.1891211* X 10~- 

Tb 

1*87.17 

1*87.17 

1*51*. 668 

325.170 

Pb 

2.302550 

1.22631*6 

.3766873 

.21061*1*0 

Pb 

.2753526 X 10“5 

.11*66537 X 10"5 

.1*826665 X 10-^ 

.3773977 x lO"' 

Hb 

1*7051.501 

5201*2.023 

61077.31*8 

79192.936 


Within the altitude range considered# T and T„ (Equation C173} 
are equal. 

Atmosphere limitations .-- The validity of the 1959 ARDC 
model IS limited to altitudes below 700 km; although the pro- 
gram is arranged to extrapolate the relationships to greater 
altitudes# if desired. Extrapolation to greater altitudes is 
accomplished by altering the Cutoff altitude. 

At an altitude greater than 2.6 x 10^ feet# no calculations 
are made, and the program sets kinematic viscosity# speed of 
sound# pressure# temperature# and density to zero. At and below 
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sea level the parameters, pressure, temperature, and density 
are set to the values below. Other terms are computed as normal. 


Pressure 

= 2116.2170 

Lb/Ft ^ 

(180) 

Temperature 

= 518.688 

®R 

(181) 

Density 

= 2.37692 X 

10“^ Slugs/Ft^ 

(182) 


At altitudes between 90 kilometers and 2.6 x 10 feet, the 
speed of sound is set to 846.50255, and kinematic viscosity 
is set to 2,3519252 x 10"^ over density. Other terms are 
computed as normal. 

The 1962 model is limited to altitudes below 295,8000 feet 
(90 kilometers) , It is suggested that zero values be returned 
above that altitude. At and below sea level, the sea level 
values should be employed. When the atmosphere constants are 
determined from the published tabulations at the base altitude, 
the calculated values at intermediate altitudes may not agree 
with the tabulated values to the number of significant figures 
in the tables. This has been allowed for in the 1959 model 
by developing coefficients with the necessary extra precision 
to give agreement between the calculated values and published 
tables at all altitudes. The values calculated by the 1962 
model are good to about four significant figures, which should 
be adequate for most purposes. 

Kinematic viscosity and speed of sound lose their physical 
significance at very high altitudes, and are not normally 
defined by model atmospheres above 90 kilometers. The constant 
values by the 1959 model option were added to provide data 
required by the aerodynamic heating routine. The aerodynamic 
heating calculation should not be used with the 1962 model 
option above 90 kilometers. The constant values of v and Vg 
in the 1959 model will give reasonable values of Mach number 
and Reynolds number for use in the aerodynamics calculations 
to altitudes somewhat above 90 kilometers, say 350,000 feet, 
above which constant aerodynamic coefficients should be used. 


Winds Aloft 

The winds-aloft subprogram provides for three separate 
methods of introducing the wind vector: as a function of 

altitude, a function of range, and a function of time. This 
facilitates the investigation of wind effects for the conven- 
tional performance studies. The wind Vector is approximated 
by a series of straight line segments for each of the methods 
mentioned above. 
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Four options are used to define the wind vector in the 
computer program. The three components of the wind vector 
in a geodetic horizon coordinate system can be specified as 
tabular listings with linear interpolations (curve reads) 
in the following options . 

Wind options (0).-- In this option the wind vector is 
zero throughout the problem. This allows the analyst the 
option of evaluating performance without the effects of wind. 

This option causes the winds-aloft computations to be bypassed. 

Wind option (1) . — In this option the components of the 
wind vector are specified as a function of time. Wind speeds 
are specified in feet per second and time in seconds. 

Wind option (2) . — The three components of the wind vector 
are introduced as a function of altitude in this option. Wind 
speed is specified in feet per second and altitude in feet. 

Wind option (3 ) . — In this option the components of the 
wind vector are introduced as a function of range. Wind speed 
is specified in feet per second and ran^e in nautical miles. 

The range utilized in this computation is the great-circle 
range . 

By staging of the wind option, it is possible to switch 
from one method <?f reading wind data to another during the 
computer run. Care must be exercised in this operation, however, 
as the switching will introduce sharp-edged gusts if there are 
sizeable differences in the wind vector from one option to 
another at the time of switching. This effect should be avoided 
except in cases where gust effects are being studied. 


Gravity 

This section presents the equations necessary for the 
introduction of the gravity components into the equations of 
motion. These components were determined by taking partial 
derivatives of the gravity potential equation. The potential 
equation adopted has been recommended for use in the Six- 
Degree-of-Freedom Flight-Path Study computer program by AFCRC. 
Constants for the potential equation were determined from 
References i8, 19 and 20. 

Spherical harmonics are normally used to define the grav- 
ity potential field of the Earth, References 18 through 2^ 

Each harmonic term in the potential is due to a deviation of 
the potential from that of a uniform sphere. In the present 
analysis the second-, third-, and fourth-order terms are con- 
sidered. The first-order term, which would account for the 
error introduced by assuming that the mass center of the Earth 
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is at the origin of the geocentric coordinate system is assumed 
to be zero. With this assumption 



where P2» P3/ and P4 are Legendre functions of geocentric lat- 
itude <pj^ expressed as 

?2 * 1 T 3 sin^ 

= 3 sin <J>L - 5 sin^ (|ij^ 

p ^ 

P 4 = 3 - 30 sin*^ 4 il + 35 sin (184) 


The gravitational acceleration along any line is the partial 
derivative of U along that line. At this point, it should be 
noted that the three mutually perpendicular directions in the 
spherical cgordinate system are identical (other than sign) to 
those in the local-geocentric-horizon coordinate system which 
is defined previously. Therefore, the acceleration in the 
direction is identical to gx„f and the acceleration in the R 
direction is identical to -g| . Or in the equation form; 


Sz 


8 


JU _ Mg. 
8R ~ ” R 




«X 


g 


1 ^ 
R ■ p2 




(-6 sin 4) cos 


k\ (185) 


(186) 


H / \ ^ / ? X 

F V W” / (3 cos ()>L - 15 sln^ (j>L cos ^^) 

/ 

I ^ I (-60 sin 4 >l cos (j)L + l4o sin3 cos ij>£) 


30 


Collecting terms: 


sz. 


8 


-ilfiL 

«2 


ex 


g 




(187) 

(1S8) 
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where 


i 


Pg = sin COB (f)^ 

» 

Pg s COS (fij^ (1-5 sin^ 

P^ = sin cos (P^ (-3 + 7 sin^ (j)^^) (189) 

Equations X13Z) ' and •(188)’ are used in the gravity subroutine 
with the following values recommended for the constants: 

Hg = 1.407698 X 10^^ ft.Vsec.^ 

» 20,925,631. ft. 

J - 1623.41 X 10-6 

K = 6.37 X 10“^ (190) 

It should be noted that these constants and equations per- 
'tain to the planet Earth; however, it is possible to use these 
same equations for any other planet. For this reason, the 
values of these constants is an input to the program so that 
the applicable constants may be inserted for the planet under 
consideration. Due to limited knowledge of the gravitational 
fields of other planets, it is probable that zero values would 
be assigned to some of the harmonic coefficients when the pro- 
gram is used for entry studies on other planets. 

The above equations are applicable to a non-rotating planet 
as the centrifugal relieving effects caused by the planet's 
rotation are included in the equations of motion. In addition, 
the effects of local anomalies must be added if it is desired 
to make a weight-to-mass conversion based on a measured weight. 
The program has the options of retaining the first, third, and 
fourth order terms. 
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AIRCRAFT CHARACTERISTICS 


• Weights 

During the 1962 F4H-1 time-to-climb record flights. Reference 21, 
a considerable effort was made to reduce the vehicle flight weight. 
Production F4H-1 and record flight vehicle empty weights were; 


Vehicle 


Empty Weight 


Production F4H-1 


29365 

- 

Record Flight Vehicle 

(3KM-15KM) 

25903 

3462 

Record Flight Vehicle 

(20KM-30KM) 

25520 

3845 


In addition to weight saving measures the record breaking vehicles 
employed non-standard General Electric J79-GE-8 engines and the record 
attempts were deliberately made during cold weather conditions to improve 
propulsion system performance characteristics for weight comparison 
purposes . 

In standard F4H-1 aircraft empty weight quoted in the 1960 weight 
statement, Reference 22, and the currently quoted empty weight of in- 
service F4-C aircraft. Reference 23, at Edwards AFB are: 

Vehicle Empty Weight 

F4H-1 (F4B) (1960 weight statement) 27345 

F4-C (Edwards AFB 1975) 32500-33500 (Net) 


For the present report it will be assumed that an average Edwards AFB F4-C 
might be made available to NASA and an empty weight (net) of 33000 lbs. 
will be utilized. In addition, Edwards AFB pilots have indicated that safety 
considerations dictate 2000 lbs. fuel on-board at approach and a require- 
ment for 600 lbs. of fuel for a descent from 40000 feet. Early calculations 
indicate that zoom climbs to high altitude will consume approximately 1600 lbs. 
of fuel. Based on these figures it is assumed that vehicle weight at zoom 
commencement should be approximately 37500 lbs. (1166 slugs). 


Propulsion 

The F4-C Phantom is powered by two J79-GE-12 engines having a 
nominal rating of 17000 lbs. static thrust at sea level per engine. Thrust 
produced by these engines and their fuel flow varies strongly with Mach 
number, altitude, and throttle setting. That is, the thrust and fuel flow 
are of the form 
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where 


T = T (M, h, N) 
W = W (M, h, N) 


(1913 

C192) 


T = Thrust, lbs. 

W = Fuel flow, lbs. /hour 
M = Mach number 
h - Altitude, feet 
N = Throttle parameter 


The F4-C engines have a wide range of throttle setting capability 
both with and without afterburners ignited. This capability is illustrated 
by Table I which presents thrust and specific fuel consumption for the 
related J79-GE-17 engine which powers the Phantom F4~E vehicle. 

Actual thrust and fuel flow capability of the J79-GE-12 engine 
which powers the F4-C is presented in tabular fashion in Table II. These 
tables are in a form acceptable to the ATOP program of References 1 and 2 
which is used to perform the trajectory optimization studies of this 
report. Format of Tables 11(a) and (b) is as follows: 

Thrust 


TTABoO = H2, 


Mj, Mj. - 




- - - 


T, 


W2 


Fuel Flow 


TTABll = H^, H^, 

Mj, M^, - - - 


Nj, Nj. 


w, 


Wi 


- w, 


N. 




M 


The data of Tables 11(a) and 11(b) was made available through the 
studies of Reference 3. 


61 


TABLE I. ’HPICAL EFFECT OF THROTTLE SETTING ON THRUST AND 

SPECIFIC FUEL CONSUMPTION OF THE PHANTOM F4-E (J79-GE-17) ENGINES 


fibs, fuel/hour/lb. thrust) 


Power 

Thrust (lbs.) 

R.P.M. 

Specific 

Fuel 

Consumption 

Max. Afterburner 

17900. 

7460 

1.965 

Military 

11870. 

7460 

.84 

Normal 

11100. 

7435 

.81 

90% Normal 

10000. 

7140 

.79 

75% Normal 

8330. 

6900 

.76 

Idle 

350. 

5000 

1130 lbs. /hour 




0-0 0 
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TABLE 11(a). THRUST OF THE J79-GE-12 ENGINE WHICH POWERS THE 

F4-C PHANTOM (cont'd). 


TTARiO 

TTARlO 

TTARIO 

ttapio 
TTARIO 
TTA»10 
T lAPlo 
TTARIO 
TTARIO 
TTARiO 
TTARIO 

TTARIO 
TTARIO 
TTARIO 
TTARIO 
TTARiO 
TTARIO 
TTARIO 
TTARIO 
TTARiO 
TTARiO 
TTARIO 
TTARIO 
TTARIO 
TTARIO 
TTARiO 
TTARIO 
T Tap 10 

ttapio 

TTARIO 
TTARIO 
TTARIO 
TTARIO 
TTARIO 
TTAPIO 
TTARIO 
TTARIO 
TTARiO 
TTARIO 
TTARIO 
TTARiO 
TTAPIO 
TTARIO 
TTARIO 
TTAPIO 
TTARIO 
TTAPIO 
TTAR 1 0 
TTARIO 
TTARIO 
TTARIO 

ttapio 
ttabio 
TAPlo 


l?aOO*., 10500',,8T30.,71R0. 
68no*,*;oiO;, 0310',, 3310', 

2500., iflpo', , ta'io'., 1100, 

Rio’., 0,0 

19R00,, lB‘J/0, , 16600,, 1«6*)0, 
1?690,, 10//U'. ,*>000*, ,7«00‘, 
7070*-, •;fllO,,<ld60,, 5420’, 

2590., iP'iO,, 1500',, 1 130, 

Bao’,, 0,0 

I93fl0'., 1*>3R0,, 176R0,, 15680, 

15550., I 1630,, 98?0',, 8080', 

7730* 6300'., OQPO,, 5760', 

2850.. 2100., 1030,, 1200. 

-p?«‘.,o,p 

0, 19900., 19510,, 17030, 

1 a 79o'., 12900., 10990,, 9060, 
POPO’-, 7120., 5520',, 62aO*, 
3220^,2920',, 1850,, 1900, 

1050. . 0.0 

0,20500., 19920., 17390, 
15180*., 13290’, , 1 1300,, 9310, 
8960*-, 7350’., 5700’,, 9380'. 

5590 ,2510., 1900,, 1950. 

1100. . 0.0 

0, 0, 19950, , 18990. 

16990*., 19350'., 12350,, 10220, 
989 0* ,5080', ,6290,, 9 850, 

3710 ,2800’., 2100,, 1610. 

1230. . 0.0 

0, 0,20300., 19630', 

178 70',, 15630'., 1369 0,, 1 1900. 
1C 960., 90 30., 7030*,, 5950, 
9190* ,3180. ,2390,, 1820', 

1390. . 0.0 
0,0,0*19800. 

18520., 16990',, 1959 0,, 12280. 
1 IPOO., 9720',, 7580'., 5890, 

959 0* ,3950’, ,2600’., 1900, 

1510. . 0.0 

0 , 0 , 0,0 

19010.. 17380* , 15550,, 13380, 

12850., 10590., 8290'., 6990'. 
9980* ,3790'., 2870’., 2170. 

1650. . 0.0 

0 , 0 , 0,0 

19200., 17 790., 16300,, 1 9 39 0, 
13P20.', 1 1910'., 8990, ,6960, 

59 00’,,9130’. , 5120*. ,2350’, 

1790. . 0.0 

0 , 0 , 0 , 0 ^ 

0 , 17900’., 16.Sno'., 19910. 

1 455o'. ,1201 o’, , 9380 , , 7360’, 
56*>0*-, 9390',, 3330,, 2500, 








TTA810 

0 » 0 # 0 » 0 


ttapi u 

0,0, 16710., 15250, 


TfAHlo 

16980*., 123*)0’.,<)700,, 7610 

1, 

TTAPlo 

5900*^ , 6560', , 3«9o; , 2620’, 


T TAP 10 

I960., 0,0 


T7AH16 

0, 0, 0, 0 


TTARIO 

0,8, 16710,, 1S370, 


TTAhlO 

150 7 o'., 12500', ,98?0’,, 7690 


TTA81U 

600 0,,6 030., 3S70*,,2680; 


TTAR16 

2010. ,0,0 


TTAPly 

0 , 0 , 0 , 0 * ^ 


TTARIu 

U,0, 16650, , 15230, 


TTARlo 

UP9o’., 12390', ,976 0'., 7630 


TTARlo 

5960* , 6620'. , ?560’. ,2680. 


TTAMo 

2020. ,0,0 


TTA^lO 

0, 0, 0, 0 


7TAR10 

0,0,0,16920. 


TTAPlo 

16510., 12060*., 9510. ,7620 


TTAPlo 

5810^, 6530'., 36/0’., 2660'. 


T tab 10 

2000, ,0,0 


TTAPlo 

0 , 0 , 0 , 0 

; J 

TTAPlo 

0 , 0 , 0 , 0 


TTAttlO 

1 5780., 1 1670’., 9060'., 7100 


TTAPlo 

5560* ,6320’. , iSiO, , 2520’. 


TTAPlo 

1950. , 0,0 


TTAPI 0 

320'.. 280’., 250’., 0 


TTAPlo 

0 , 0 , 0 , 0 


TTAPlo 

0 , 0 , 0 , 0 


TTABlO 

0, 0, 0, 0 


TTAPlo 


i 

TTAPlo 

260. ,220., 200. ,0 

) 

t 

TTAPlo 

0 , 0 , 0 » 0 


TTABlO 

0 , 0 , 0 » 0 


TTABlO 

0 , 0 , 0 , 0 


TTAPlo 

■ 0,0,0 ■ 


TTAPlo 

200., IflO*., 16O’.,0 


TTAPlo 

0 , 0 , 0 , 0 


TTAPlo 

0,0, 0,0 


TTAPlo 

0 , 0 , 0 , 0 

> 

TTAPlo 

0*0/0 


TTAPlo 

•60. ,-30., -10. ,0 


TTAPlo 

10*., 30’., 1 10. ,200, 


TTAPlo 

230., 300'., 350. ,350', 

j 

TTAPlo 

0 , 0 , 0 , 0 


TTAPlo 

0,0,0 

o 

TTAPlo 

•310l,-260, ,-210.,-l65. 


TTAPlo 

-130.,-75.,1C.,105. 


TTAPlo 

120’., 190*., 255. ,290. 


TTAPlo 

0 , 0 , 0 , 0 

• 

TTAPlo 

0,0,0 


TTAPlo 

-610, -510, ,-660, ,-355. 

: 

T T P i 0 

-295.,.220..-t25,,-lo’, 

i 

TTAPlo 

10*., 70’., 155. ,220', 


T 


•\27 
Oil 
035 
030 
00 3 
006 
050 
050 
058 
062 
065 
060 
073 
077 
Ofll 
08 0 
088 
002 
006 
500 
503 
50 7 
511 
515 
510 
522 
*■>26 
530 
5 30 
536 
50 1 
505 

500 
55 3 
557 
560 
560 
568 
5 72 
5 76 
570 
583 
587 

501 
505 
508 
602 
606 
610 
610 
617 
621 
625 
620 


Hcl.OSO ‘ 

i 

8»2,ino • 
8«2’,20C 
M»2,3no 
MB?, 000 
IOLC»>*aO 

i 

*<»,05o 

I 

*^s,l00 

8»,250 

Ms,o0o 

M*,55n 


V 




u 


u 


o 


o 


TTARIu 

TTAHIO 

TTAPIO 

TTARlCi 

TTARIu 

TTARIu 

TTARIu 

ttariq 

TTARIO 

TTARlO 

TTARIu 

TTARlO 

TTARIu 

TTARIu 

TTARlO 

TTARlO 

TTARIu 

TTARIU 

TTARlO 

TTARlO 

TTARIu 

TTARlO 

TTARIu 

TTARlO 

TTARIu 

TTARlO 

TTARlO 

TTARIu 

TTARIu 

TTARIu 

TTARlO 

TTARIu 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARIu 

TTARlO 

TTARlO 

TTARlO 

TTAPIO 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARlO 

TTARIu 

TTARlO 

TTARlO 


T»»0,O, 

• 100S,»"860,f««730,»»620t 
-5io',,-'SRo.,-2/*o;,-tao; 
-100. ,-00, ,60., ISO*. 

0 , 0 , 0 , 0 


0 , 0,0 

-1380*. ,-l 1 •50., -1 000,, -860, 
-6^0. ,-52^, ,-370. ,-220, 
-lRO.,-100,.0, lOO, 

0 1 0 , 0 , 0 


0 

0 , 0,0 


-855;, 

-655', , -SUO’, , - 

•280., 

-190'., -60. ,50 

0,0,0, 

0 


0,0,0 



-2050; 

-1095. 

,-2010* 
, —880 , 1 

,,-1600 

,-695*., 

-005'., 

-310.,. 

•160,,- 


0 , w , w , W 

0 , 0,0 

•230o' ,-2llO.,-170‘5,,-1010, 
•l|RO.,-OtO,, -725’, , -505, 
.a 75 ’.,. 300 '.,-t 80 ,,- 5 o’. 

), 0,0,0 


0 , r M « 

0 , 0,0 

R0 0’.,-2i'30,,-l9l0.,-l560, 
1250. ,-1010. ,-810’., -600’, 

550’.,-005;,-230,,-90. 

I . ft . n - 0 


0 


-550., 
0 , 0 , 0,0 
o,o»o . 

,-790., 

1 « 20 ’ 
055*. , 

0, 0, 0 , 0 

0 , 0,0 


'90’.,-?lflo'.,-l755. 

;o’., -1100’., -910, ,-700. 
i*.,-o95‘.,-300’,,-ia5. 

A ft 


0 , 0,0 

0, -100. ,-2050., -1810, 
-l«60’. ,-1160’,, -900’. , -725, 
-69O*.,-520’.,-3l0‘.,-l55, 

0 , 0 , 0 , 0 


0 , 0,0 

0,0,0.-1960. 

-1590., -1260’, ,-1015., -795, 
-755'. ,-600., -375*., -210, 
0 , 0 , 0,0 
0 , 0,0 

0,0, 3000’., 600'. 

-l55o’. ,-1390’, ,-l09u,,-855, 
-810. ,— 670.,— fl50.,— 275, 

0, 0, u, 0 
0 , 0,0 

0,0,0,2650’. 

350. ,-1200, ,-1125., -I 
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63} 
6iU 
6a 0 
6a/i 
6'i8 
652 
655 
659 
663 
667 
671 
6 fn 
678 
682 
686 
690 
093 
697 
701 
705 
709 
712 
716 
720 
720 
728 
731 
735 
739 
7a3 
7‘17 
750 
75a 
758 
762 
766 
769 

773 

777 

781 

785 

708 

792 

796 

800 

8oa 

807 

811 

815 

810 

823 

826 

830 

859 


MB, 700 


M«,800 


Ma,900 


Moi ,050 


M»l ,0/7 


M»1,150 


MBl‘,2a7 


hbI ,272 


MBJ,360 


M»1 ,a/o 


MBl ,550 



! 


TABLE 11(a). THRUST OF THE J79-GE- 12 ENGINE WHICH POWERS THE 

F4-C PHANTOM (cont'd). 


TTAPlO 

0, 0, 0, 0 

TTAPIO 

0,0,0 

TTARlO 

0 , 0 , 0 , 0 

TTARlO 

?a00* 525, ,-11/5, ,-925’, 

TTARlO 

-R70.,-A90’.,-595.,-395. 

TTARlO 

0 , 0 , 0 , 0 

TTARlO 

0,0,0 

TTARlO 

0 , 0 , p , 0 

TTARlO 

12R0,,2/|50-,660,,«»955, 

TTARlO 

-90 0.,-780.,-650’. ,-«70’, 

TTAHlo 

0 , 0 , 0 , 0 

TTARlO 

0,0,0 

TTARlO 

6, 0,0,0 

TTARlO 

0,1 2350’., 510*. 

TTARlO 

RO’., 130., 75. ,60, 

TTARlO 

0 ,• 0 , 0 , 0 

TTARlO 

0,0,0 

TTARlO 

0 , c , 0 , 0 

TTARlO 

0,0,1300*. ,2^100'. 

TTARlO 

2300*., 1950'., 1500’,, 1200’, 

TTARlO 

0, 0, 0, 0 

TTARlO 

0,0,0 

TTARlO 

0, 0, 0, 0 

TTARlO 

0,0,210*., 1 170, 

TTARlO 

1275'., 1030-., Rio, ,64J5*. 

TTABI 0 

0,0, 0,0 

TTARlO 

0,0,0 

TTARlO 

0, 0, 0, 0 

TTAPIO 

0,0.-700* , 390', 

TTARlO 

575,,d70., 370.,2fl0'. 

TTARlO 

0, 0,0,0 

TTARlO 

0,0,0 

TTARlO 

0, 0,0,0 

TTARlO 

0, n, 0, -060’, 

TTARlO 

-?30., -180,, -1^0, ,-125, 

TTARlO 

0, 0, C, 0 

TTARlO 

0,0,0 

TTARlO 

0 , 0 , 0 , 0 

TTARlO 

0 , 0 , 0 ^ 0 

TTARlO 

-l240,,-l000,,-300.,— 640 

TTARlO 

0 , 0 , 0 , 0 

TTARlO 

0,0,0 




mbj.650 


MSI ’,750 " 




►’Cl ,950 


MS?’, loo 


Mc ?,?00 


Ms ?,«00 


o 
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TABLE 11(b). FUEL FLOW OF THE J79-GE-12 ENGINE WHICH 
?S 

2 THE F4-C PHANTOM. 

o,*’»0uo., lonon,, isooo’. 

r^'OOo' ,2*3000', ,30000,. 5500 0, 
iOObQ ,00000, , 05000, ,50000, 

55000, ,6000 0, ,650 00,^7 0000', 

75000., 75010, , 120000, 

0,*,050,'. 100,'.250 
>00,*.550,',700, ,fl00 

^ n IS n * * rn . % ' t Kn 


2fl«00;,20600, ,20700,,0 
0 , 0 , 0 , 0 
0,0, 0, 0 
0 , 0 , 0 , 0 
0,0,0 

30 75o'., 2630 0',, 220 00,. 0 
0, 0, 0, 0 
0 , 0 , 0 , 0 
0,0, 0,0 
0,0,0 

31650*., 27050',, 22650, ,0 
0, 0, 0, 0 
0, 0, 0, 0 
0 , 0 , 0 , 0 
0,O,0 

330P o',, 29600*,, 20 120,, 20050, 

16750.. 0.0.0 

0,0, 0,0 

0, 0, 0, 0 
0,0,0 

35300*,, 30390'., 25730,, 2 1500, 

17950., 10700., 12190,, <>950 , 
0590*. ,0, 0,0 

o,r,o,o 

0,),0. 

37500* ,32600' 27950,, 23050, 

10100.. 15050.. 13080.. 10790, 

10310.. 6750. .6,)ao'.,0 
0 , 0 , 0 , 0 

0,0,0 

4 05 70*,, 353201, 30550,, 25060', 
21 600 ;, 17750 , H55o., 11070, 

11390. . 0630. .7700,., 6030. 
4690*., 0,0,0 

0,0,0 

43290*,, 37050', ,32550., 270?O‘, 

23420., 1 '>350,, 15920,, 1206O', 
1231 0., 10390,, 9300’,, 6570', 
5090*., 3970',, 3060',, 2360, 
0,0,0 

46550*, , 40260', , 34 75o, , 290?O , 

25380.. 2 11 30.. 17350.. 14050, 

13380., 1 1280'., 9110'., 7180', 

55 70*., 43 10'., 3350,, 2590', 
0,0,0 

4 765o', , 45450'. , 3890O, , 33420 . 

29550., 241 0 0'. ,20 020,, 16250, 

15550.. 13040.. 10490.. 9340, 


o 


o 


TABLE II (b). FUEL FLOW OF THE J79-GE-12 ENGINE WHICH POl^fERS 

THE F4-C PHANTOM (cont'd). 


\o 


i O 


o 


o 


TTAPIi* 

29100' ,20050*. ,20550,, 16670', 

261 

T TAPI 1 

1601 0., 13020’. , 10 / 60 ,, 8600, 

265 

TTARl 1 

6090* ,5170'. ,« 000'., 3090'. 

269 

TTARU 

2380., 0,0 

253 

TTABU 

a7 170,, '16650’, ,«2220,, 36380’, 

256 

TTAPlj 

30/50,, 26260., 22000., 17950, 

260 

TTA«l 1 

17270., lOOSO., 1 15/0. ,9270. 

266 

TTAPl 1 

7230* ,560 0*. ,6 330., 336 O’, 

268 

ttahu 

2570. ,0,0 

272 

TTAPl 1 

0, 6 /JOn., 65850., 39530, 

2/5 

TTABI I 

3335o’.,2R650'. ,2600 0,, 19770, 

2/9 

TTAPII 

19020., 15900'., 12/10,, 10210. 

283 

TTAHll 

8010*,, t250‘.,«fl20*., 3720’. 

207 

TTAPII 

. 2880. ,0,0 

291 

TTAPII 

0,67330. ,65980’., 60320'. 

296 

ttapu 

36 1 Oo', , 29120’. , 26570, , 20300’, 

29R 

TTA^H 

19520., 16310, , 13020,, 10680. 

302 

TTaHI 1 

8270* ,6620’., 6970, ,3800'. 

306 

T TAB 11 

2970. ,0,0 

310 

TTaPIJ 

(I, 0,66370. ,62680’. 

313 

TTABll 

37! 50*., 3 1680', ,26550,, 22030’, 

317 

TTAfll 1 

21 l8o., 1 7690’. , 16100,, 11320. 

321 

TTaPI 1 

9020,, 7060’., 56/0, ,6250, 

325 

TTAWI 1 

3290. ,0,0 

329 

TTAPII 

0,0, 66670, ,65000’, 

332 

TTABi j 

6 0250*,, 36750'., 29200,, 26350, 

336 

TTAPII 

25350,, 195 SO,, 15520., 12630, 

360 

ttarii 

100 00. , 7900, , 6120’, ,6/70, 

316 

TTAPIi 

3700*., 0,0 

560 

TTARl I 

0,0,0.65620. 

351 

TTABfi 1 

6190 0,, 36650',, 3 1330,, 26 150'. 

355 

TTARl 1 

25070,, 20960., 16630., 13310, 

359 

TTARII 

1 0680., 85/0., 6630’., 5180, 

363 

TTARl 1 

6000*., 0,0 

567 

TTAPU 

0, 0, 0. 0 

370 

ttarii 

63600,, 38500*., 339/0,, 28650', 

3 76 

TTAPII 

27360,, 227 SO., 13090., 16620. 

378 

TTAPII 

1 1550. ,9330. ,7330*., 5700. 

382 

TTAPIi 

6620*. ,0,0 

386 

TTARl i 

0 , 0 , 0 , 0 

589 

TTABll 

66 3 30,, 39850’, ,35/70,, 3070 0, 

393 

TTARl i 

29650,, 26600., 19560,, 15580, 

397 

TTARII 

12650., ICO SO'., 8030’., 6270, 

601 

TTABI i 

6870*., 0,0 

605 

TTARl i 

0,0, 0,0^ 

608 

TTARII 

0,60670., 36670, i32300. 

612 

TTAPII 

3131 0*,, 26050'., 20670,, 16660’, 

6 1 6 

ttarii 

1 3100., 10520'. ,85 10'., 6710, 

620 

ttarii 

5210’., 0,0 

626 

TTARl 1 

0, 0, 0, 0 

627 

ttarii 

0,0, 37600., 33670’. 69 

631 

TTARl 1 

32620*,, 27190’,, 2 1560,, 17120’. 

635 

TTABll 

13660^, 1 0930’., 8860’., 71 00', 

639 

TTABll 

5560*.. 0,0 

663 








.“SOO 


HBl ,«/o 


M*1 ,?)50 I 


HsitbOO . I 


«»1 ,750 


,8d0 


mbi ,950 


j 
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I TABLE 11(b). FUEL FLOW OF THE J79-GE- 12 ENGINE WHICH POWERS 
I THE F4-C PHANTOM (cont'd). 

o 


u 




i ) 


|o 


iG 


TTABl i 

Of Of Of u 

TTAPI 1 

0,0, 38330, ,3aa5C, 

TTAPI 1 

33680^,28080^,22220, 

TTARl 1 

14100., 1 1260, ,0110, ( 

TTARll 

5820*., 0,0 

T TAPI ^ 

0,0, 0,0 

TTARll 

0,0, 38000, ,34800', 

TTARll 

33O30* ,28300' ,22440, 

TTARll 

14200., 11380. ,0190,, 

TTARll 

5040*., 0,0 

TTARll 

0 , 0 , 0 , 0 

TTARll 

0,0,0,30050, 

TTARll 

33820-, 28100’,, 22380 

TTARll 

1417U., 1 1 370’., 0160'., 

TTARll 

5000*., 0,0 

TTARll 

0 , 0 , 0 , 0 

TTARll 

0, 0, 0, 0 

TTARl 1 

33140-, 27640’., 21070 

TTARll 

I30iu., 1 1 160. ,0020', 

TTAPI 1 

5060*., 0,0 

TTARll 

1275*., 1085’., 0l5,,0 

TTARll 

0, 0, 0, 0 

TTABl 1 

0 , 0 , 0 , 0 

TTARll 

0, 0, U, 0 

TTARll 

0,0,0 

TTARl 1 

t 105., 1005’., 860,, 0 

TTAPI i 

0 , 0 , 0 , 0 

TTARll 

0, 0, 0, 0 

TTARll 

0 , 0 , 0 , 0 

TTARll 

0,0,0 

TTARll 

1140*., 070., 835’,, 0 

TTABl 1 

0 , 0 , 0 , 0 

TTARll 

0 , 0 , 0 , 0 

TTAPI 1 

0,0, 0,0 

TTAR J 1 

0,0,0 

TTARll 

1045*., 000., 765'., 655 

TTARl 1 

550'^, 500*-, 500., 500, 

TTAPI 1 

500, ,500, ,500. ,500. 

TTARl 1 

0, 0, 0, 0 

TTARl 1 

0,0,0 

TTAPI 1 

055, ,810*, 700., 60C, 

TTAPll 

500,,500,,500,,500: 

TTARll 

500. ,500, ,500. ,500. 

TTAPI 1 

0, 0, 0, 0 

TTARl 1 

0,0,0 

# A ^ , ■ A 

TTARll 

80S ,685,, 58S., 500, 

TTARll 

500 ,500, ,500. ,500, 

TTARll 

500. ,500, ,500. ,500'. 


TTAPI 1 

ttapi 1 
Tr»oi i 
TTAPI i 
TTAPU 


7370’, 


7^30’. 


, 1 

7320, 


U| 0, 0 


r.7S',.«5OO’-,5OO.,SO0 
•500,,*500,,S00 
500 
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/|<I6 
050 
05^1 
058 
062 
065 
060 
0 73 
0 77 
081 
080 
088 
002 
00b 

500 
50 3 

507 
5U 
515 
510 
522 
526 
530 
530 
538 

501 
50 5 

500 
553 
557 
560 
560 
568 
572 
576 
5 70 
583 
587 

501 
505 

508 
602 
606 
610 
610 
617 
621 
625 
620 
633 
636 
640 
640 
646 


mb?, 1 OO 


M»2,?00 

i 

I 

mb?',3o0 ‘ 


MB ?,«00 

\ 

I 

i 

i 

l0UC,«*0i 


mb,050 


mb, 100 


mb, 250 I 

I 

I 

mb,«00 


Ms,55o 


MS, 700 



TABLE 11(b). FUEL FLOW OF THE J7&-GE-12 ENGINE WHICH POWERS j 

THE F4-C PHANTOM (cont'd). ! 


0 , 0,0 

500’ •;rj‘,,*‘»oo.,5oo, 

500 ,500.,500.,500j, 

500. . 500. .500. .500. 

0 , 0 , 0,0 

0 . 0,0 

1305.. 1000.. 005. .6.5. 
565'^, 500’^, 500^, 500', 

500. . 500. .500. .500. 

0 , 0 , 0 , 0 

0 , 0,0 

4 *>AA* 44CC* 


~ w w - 

0 , 0 , 0 , 0 
0 , 0,0 

2000’., 1025. ,1105. ,985. 
AfiO' ,650',,5CO. ,500*. 

500.. 500..5CO.,500, 

U, 0, 0, 0 

0 , 0,0 

«375., 1555., 1205., 1070. 
87C-,715;,550.,500, 

500. . 500. .500. .500. 

0 , 0 , 0 , 0 

0 , 0,0 

O,d030'., n«5'., 1200', 

990’ ,805' ,640. ,500’ 

500. ,500. ,500. ,500. 

0 , 0 , 0 , 0 

0 , 0,0 , 

0,d900., 1/20. ,1235. 
1015*.,835.,6S5.,52 o, 
500*., 500'., 500. ,500. 

0 , 0 , 0,0 
0 , 0,0 . 

0,8500. ,0050., 1000. 
1105*., 920, ,735. ,585. 
545'., 500’., SOO'., 500, 

0, 0, 0,0 
0,0,0 

0,0,8500*. ,0600, 

1 0 00*. , 1055‘. , 805, ,6/5*. 
630*., 505*., 500. ,500’, 

0 , 0 , 0 , 0 
0 , 0,0 

0,0,0,7000. ^ 

i860*., 1390’. ,945, ,755. 
700’., 575'., 500’., 500, 

0 , 0 , 0 , 0 
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TABLE 11(b). FUEL FLOW OF THE J79-GE- 12 ENGINE WHICH POWERS 1 


THE F4-C PHANTOM (cont'd). 


ttapi 1 
TTAP-U 
ttapu 
tta«u 
ttapi 1 

TTAP 1 i 

ttapi 1 
ttapi i 

TTaP I \ 
TTAPU 

ttaru 
T r ah u 
ttah u 

TTAH I 1 

ttabu 

TTAHI i 

ttapu 
tta«u 
ttapi 1 
ttapi 1 

T T A p 1 i 

ttapu 
ttapi I 
ttapi i 

TTAPU 

ttapu 

TTAR) \ 

ttapu 

T TA«1 1 

ttapu 
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Aerodynamics 


Aerodynamics for the W-C are functions of Mach number, altitude, 
and angle-of-attack. For three-degree-of-freedom trajectory calculations 
lift and drag coefficient variations are required. Lift coefficient is 
a direct table lookup of the form 

(a, M, h) 

Drag coefficient has two components 

Cjj - (a, M, h) + (M) 

Actual values of lift coefficient and drag coefficients employed for the 
F4-C trajectory calculations are given in Tables III (a) to lll(c). 

Tabular format of these tables is similar to that employed in the pro- 
pulsion section. Source of the F4-C data is the Air Force Flight Dynamics 
Laboratory through the study of Reference 9. 
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TABLE Ill(b). DRAG COEFFICIENT vs. Cj^ and MACH , F4-C AIRCRAFT 
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TABLE III(c). ADDITIONAL DRAG COEFFICIENT vs. MACH NUMBER 

F4-C AIRCRAFT 
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DATA CALIBRATION BY 
LEVEL ACCELERATIONS 


Trajectory Calculations 

A series of level acceleration trajectory computations were 
performed to delineate maximum Mach number capabi^ties using the 
aircraft propulsive and aerodynamic characteristics of the previous 
section. Trajectories were simulated by using the commanded flight 
path angle ( 7 ) option of the Reference land 2 program. That is, the 
vehicle flies an angle-of-attack schedule which attempts to maintain a 
specified flight path angle history, in this case 

7 = T(t) - 0 (195) 

Table IV(a) presents a typical acceleration starting at 40000 feet with 
an initial velocity of 1000 feet per second. Vehicle initial mass is 
1250 slugs (40217 lbs. weight). Maximum Mach number achieved is M = 2.06 
at a time T = 375 seconds and range R = 99.02 nautical miles. Theoretical 
maximum altitude capability at the end point based on energy conversion is 
E = 102000 feet. Table IV(b) presents the terminal states achieved in a 
series of level accelerations starting at initial altitudes varying from 
20000 to 50000 feet. 


Some points should be ma4e regarding these trajectory calculations. 
First , the maximum Mach number achievable is a function of the vehicle 
weight. As the weight diminishes vehicle lift coefficient for level flight 
diminishes and hence a reduction Occurs in drag-due-to-lift . As the drag 
diminishes the vehicle speed can increase until thrust again equals drag. 
Maximum speed therefore increases slightly with reduced vehicle weight. 
Second , the flight path control used attempts to steer a given flight path 
angle by an iterative numerical process. Over the period of time involved 
(approximately 400 seconds) slight errors in flight path angle cause a 
cumulative altitude error. Thus initial and final altitudes differ to 
some extent. The final flight path angles of these simulations have little 
effect on maximum Mach number. Third , maximum trajectory simulation time 
was 400 seconds, _ At the highest altitude this is insufficient time to 
develop maximum velocities. Thus the maximum attainable Mach number at 
altitudes in excess of 40000 feet, and the corresponding theoretical 
altitudes are undefined by the simulations. Subsequent simulations revealed 
a maximum Mach number capability of 2.0 at 50000 feet. - This results in a 
theoretical altitude capability of 107000 feet. 

Figure 20 illustrates the level acceleration trajectories in the Mach 
altitude plane. Time checks are displayed at 50 second intervals. At the 
lower altitudes the closeness of the time checks at trajectory termination 
is an indication that the limiting Mach number for given vehicle weight 
i*: being approached. At the higher altitudes (H > 40000 feet) very slow 
acceleration capability is encountered. However, as the vehicle Mach 
number increases the increased spacing of the 50 second time checks 
indicates that the vehicle is travelling below its limiting Mach number. 


TABLE IV (a) 

EXPLORATORY F^»-C ACCELERATION 
FROM V = 1000 FT/SEC @ 40000' 


T 

M 

H (FT) 

_T° 

W (LBS) 

R (NM) 

E (FT) 

0 

1.03 

40000. 

0 

40217. 

0 

55353. 

50 

1.21 

39889. 

1.2 

39825. 

8.93 

61127. 

100 

1.40 

39908. 

1.4 

39355. 

19.30 

68307. 

150 

1.60 

40121. 

1.6 

38797. 

31.20 

77165. 

200 

1.79 

40348. 

1.8 

38148. 

44.65 

86892. 

250 

1.94 

40946. 

0.1 

37430. 

59.45 

^ 94965. 

300 

2.01 

40583. 

0.1 

36687. 

75.09 

99418. 

350 

2.05 

40628. 

0.02 

35938. 

91.11 

101421. 

375 

2.06 

40637. 

0.01 

35569. 

99.02 

101951. 


TABLE IV (b) 

TERMINAL STATES F4-C ACCELERATIONS 


RUN 


T 







I 

40000. 

375. 

2.06 

40637. 

.01 

35569. 

99.02 

101951. 

2 

40000. 

400. 

2.06 

41320. 

.02 

35366. 

105.77 

102703. 

3 

35000. 

400. 

2.06 

39257. 

.24 

34163. 

112.06 

100935. 

4 

35000. 

400. 

2.06 

39259. 

.24 

34163. 

112.08 

100936. 

S 

30000. 

400. 

2.03 

38417, 

.78 

33462. 

111.20 

98196. 

6 

25000. 

350. 

1.80 

26610. 

.26 

33057. 

94.72 

77742. 

7 

20000. 

250. 

1.47 

14048 

-1.13 

34552. 

61.54 

52006. 

8 

45000. 

Voo. 

1.91 

45778. 

.20 

36715 

93.04 

98881 

9 

50000. 

400. 

1.35 

49280. 

.12 

38062. 

73.30 

754P0. 




Maximum Mach Number Correlation 

Level acceleration flight path calculations have been performed to 
establish F4-C maximum Mach number capability in the altitude range of 
prime interest for zoom commencement. Based on the results of Reference 25 
this altitude is about 40000 feet. These calculations also serve to 
calibrate vehicle aerodynamic and propulsive characteristics employed. 

Thus in Reference 10 discussions with pilots who fly IF4-C aircraft have 
indicated a maximum Mach number capability of 2.10 on occasion. The 
present maximum Mach number of 2.06 is in close agreement with this figure. 

It appears therefore that the aerodynamic and propulsive data of the previous 
section provides a good but slightly conservative model of the F4-C 
aircraft. 

It may also be noted that Reference 26 provides a predicted Mach- 
altitude steady state flight envelope for the F4-C aircraft. This flight 
envelope is reproduced in Figure 21 for reference purposes. Again it may 
be seen that a maximum Mach number capability of 2.10 is predicted in the 
vicinity of 40000 feet. However, the vehicle weight employed in these 
calculations is not noted in Reference 26. The zoom weight of Reference 26 
is 33238 pounds, somewhat below the value of 37500 lbs. employed in the 
present calculations. 
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ENERGY MANEUVERABILITY METHOD 


In the early fifties Rutowski, Reference 27 and Lush, Reference 28, 
independently proposed an energy maneuverability method for certain 
trajectory optimization problems. Their underlying assumption was that 
kinetic (velocity) and potential (altitude) energies are readily exchange- 
able and that the major factor in climbing to a flight condition (V, h) 
is the build-up of a specific energy corresponding to the flight point. 
Specific energy is defined as the vehicle energy content per unit weight of 
vehicle. That is 


E = mg (h + V'‘/2g) 


Eg = h + VV2g 


(196) 


(197) 


where 


E = Vehicle energy 


Eg = Specific energy 


Further, accepting the hypothesis that energy build up is the dominant 
factor in achieving a given flight condition, local optimization at each 
poiiit along a trajectory leads to the conclusion that the time derivative 
of specific energy, dE /dt, should be maximized at each energy level. This 
in turn implies that tEe vehicle should follow the loci of the tangency 
points between the contours 

2 

Eg = h + V /2g = constant (198) 


Es = dEg/dt = V(Tcosa - D)/mg = constant 


(199) 


Asssoming the thrust is parallel to the velocity, the work done in time At is 


AE = (T - D) VAt 


( 200 ) 


= i limit [flj = CT - 


( 201 ) 


This last expression may also be recognized as the first-order expression 
for Steady-state rate-of-climb. Assuming equilibrium flight and small 
angles of attack 

(T - D) = W sin 7 (202) 






and it is assumed that an aircraft flown at constant energy requires 
negligible time for state changes. 



The accuracy of the energy method has been examined in some detail 
in Reference 26. Both time-to-climb and the vehicle's ability to exchange 
kinetic and potential energies were considered. |For maneuvers near the 
limit of a vehicle's zoom altitude capability time estimates by the energy 
method are inappropriate since a major part of the flight path involves a 
near constant -energy arc. Again the assumption of perfect ability to 
exchange kinetic for potential energy is seriously in error. To completely 
convert a vehicle's kinetic energy to potential energy requires vertical 
flight at the apogee. To achieve vertical flight from the higher energy 
states of a high performance aircraft involves a pullup from almost hori- 
zontal flight. This pullup creates large drag due to lift increments 
which dissipate energy. Maximum altitude capability is thus a sensitive 
trade between drag losses which dissipate energy,; engine thrust, and energy 
producing capability near the aircraft operating limits, and the tendency 
to seek a high pullup angle to maximize energy exchange. In such a situation 
more exact methods of trajectory optimization based on the variational 
calculus must be employed. All optimum trajectories in this report are 
obtained through the ATOP program of References 1 and 2. Theoretical 
altitude capability based on energy content at zoom commencement will be 
employed only as a reference measure. This idealized altitude capability 
can be attained only if enough thrust is available to balance the drag 
forces during and after the high-g pullup. 
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PROGRAM VERIFICATION 


In 1954, Rutowski published the energy-climb optimization procedure. 
Reference 27, which was based on a localized maximization of the aircraft 
energy buildup in the Mach-altitude plane. This was the first of several 
efforts which were developed in the 1960's and which used variational 
calculus methods for solving aircraft performance-optimization problems. 

The: present program is based on the use of a variational steepest-descent 
algorithm and the formulation includes stage points and branched-trajectory 
options. The program is both versatile and computationally efficient, 
particularly when used with extremely complex high-ordered system equations. 

The success of the variational steepest-descent method in the solution 
of aircraft performance optimization problems is evident from the strong 
support given to this technique by leading government research centers. The 
reason for this support is clear when the performance gains obtained by this 
method are examined. In Reference 25 are presented results related to 
several F-4B time-to-climb flights. Of particular interest are Figures 2 
and 3 of this reference, which show paths in the Mach-altitude plane for a 
number, of record flights. The actual flight paths flown are compared with 
the paths predicted by the programs of References 1 and 2, and it is found 
that a considerable increase in performance is obtained by attempting to 
follow the predicted optimum. In a typical case, a 23% improvement in 
performance over that of the flight handbook is obtained by both theory and 
actual flight, even when the path is not stringently followed. These results 
are included as Appendix A of the present report for reference purposes. 

It is noted that this gain is obtained without vehicle modification, 
but simply by controlling the aircraft in a near optimal manner. It should 
also be observed that the performance sensitivity to departures from the 
optimal flight schedule is not very large, and that the derived optimal 
path is not difficult to implement. 

Multiple Extremals 

In certain applications, it may happen that multiple extremals exist 
between the initial and final conditions, such that small perturbations of 
each extremal produce a decrease in performance. During the present study 
two distincts sets of extremal trajectories were found from the initial 
conditions selected in Figure 21. Since the time-of-flight was free, and only 
the final altitude was to be maximized, it is apparent that this maximum 
altitude should be independent of the aircraft initial condition. Nevertheless 
from both initial conditions, the F4-C immediately dives to lower altitudes, 
and this is followed by an abrupt pull-up to the near constant-energy zoom 
maneuver. The final altitude reached from the higher initial energy state 
(point II), however, was superior to that achieved on the path beginning 
at the lower energy state (point I). Since a trajectory could physically 
connect points I and II, and since the final altitude achieved from point II 
is greater than that achieved from point I, it is concluded that two extremal 
trajectories must emanate from point I. The optimization process has con- 
verged to the inferior extremal, (which in this ease is a shorter flight) 
as illustrated in Figure 22, The miiltiple extremal problem is eliminated by 
beginning the computation at the highest energy level, for then the subaxc 
I-II no longer exists. 
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NUMERICAL RESULTS 

The optimization program has been applied to the basic F4-C aircraft, 
as described by its inertial, thrust and aerodynamic characteristics. The 
parameter of principal interest was the terminal dynamic pressure, since the 
low density and low speeds at maximum altitude could make the aircraft 
difficult to control at this time. Additional results relate to the effects 
of increased thrust, horizontal winds, and decreased gross weight of the 
aircraft . 

Terminal Dynamic Pressure Limits on Altitude Capability 

The computations were performed with the F4-G initially at 50000 ft 
altitude, and at M = 2.06. The initial mass was 1166 slugs, and the terminal 
dynamic; pressure was specified at 10, 20 or 40 psf. Results are shown in 
Figure 23, and it is seen that the maximvun altitude varies from 88191 
(qf =8.7 psf) to 82780 (qf = 39.9 psf). 

The question of specifying an adequate minimum dynamic pressure is one 
which must be answered by use of detailed simulations using both transla- 
tional and rotational degrees of freedom. These simulations should include 
details as to the flight control system, and the aerodynamic stability 
deviations at large angles of attack and large sideslip angles. For 
reference purposes, it may be noted that recent F-15 high altitude time-to- 
climb record flight went safely over the top with a dynamic pressure less 
than 10 psf, Reference 29. It is doubtful that routine sampling of the 
upper atmosphere would be performed at such low dynamic pressures however. 

Effect of Increased Thrust on Maximum Altitude 

When the available thrust is increased by 10%, the maximum altitude 
increases, as shown in Figure 24. According to data provided by General 
Electric to McDonnell Aircraft, Reference 30, a 2% increase in engine rpm 
at M = 2.0 generates a thrust increase of the F4-C of 10%. This is increased 
to 12% at M = 2.2, and is decreased to 6% at M = 1.9. Since the Mach number 
on the nominal flight exceeds M = 2.0 for about 75% of the powered portion 
of the zoom climb maneuver, an average value of thrust increase equal to 
10% has been assianed. As shown in Figure 24, an increment of 2000 to 4000 ft 
altitude is obtained from this thrust increase, depending on the final 
dynamic pressure specified. 

Effect of Tail Winds on Maximum Altitude 

The Mach number attainable by the F4-C at a particular altitude is 
independent of the winds, but the kinetic energy is not. That is, the 
inertial velocity of the aircraft is increased when a tail wind is acting 
on the aircraft, and this additional inertial velocity can be converted to 
an increment in final altitude. The wind profile given in Figure 25 was 
taken as representative, and it is noted that the initial extra velocity, 
at 50000 ft altitude is only 75 ft/sec. 
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FIGURE 24. EFFECT OF EXTRA THRUST ON MAXIMUM ALTITUDE 
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The effects of this tail wind profile are shown in Figure! 26, for 
terminal dynamic pressures of 10, 20 and 40 Ib/ft^, The 3 curves indicate 
that the wind and the extra thrust generate approximately equal altitude 
increments. It is a matter of interest that the minimum altitudes reached 
during the zoom-dive vary from about 33000 ft (qf = 10 Ib/ft^) to 43000 ft 
(qf = 40 Ib/ft^), and the peak wind velocities were assumed to occur 
between 30000 and 40000 ft. 

Effect of Reduced Initial Mass on Maximum Altitude 

When the aircraft mass is reduced due to the burning of fuel, greater 
accelerations result from a given thrust, and this has a modest effect on 
the altitude capabilities. The initial weight was reduced from 37500 lb to 
32170 lb. The optimization procedure was then carried out at the initial 
Mach number of 2.10, because at this speed the drag is in equilibrium with 
the 10% additional thrust. 

As shown in Figure 27 for a final dynamic pressure of 20 psf, the 
increase in maximum altitude due to this mass reduction is 1750 ft. 


Transients in Selected Trajectory Variables 

The transient variations bf Mach number, angle of attack, altitude, 
etc., during a typical optimal flight are shown here in graphical form, to 
indicate the expected values attained by these variables. The flight chosen 
for this example is that for which the altitude reached is 92.6kft and the 
minimum dynamic pressure is 20 psf. This example includes the effects of 
winds and the additional thrust while the gross weight is at its nominal 
value. Asi shown in Figure 28, the dynamic pressure and the normal acceler- 
ation reach peak values during the zoom-maneuver, when the altitude is near 
its minimum. The angle of attack varies more smoothly, and has a peak value 
approximately mid-way between the minimum and maximum altitudes, when the 
flight path angle is also near its peak value of 45°. 

These results indicate that control variations during such a flight 
need not be abrupt or extreme, and that normal accelerations can be kept 
below 4-g'!s. In fact, the output data for the 12 cases given in Table V 
shows that the largest normal acceleration encountered is 3.96 g’s, which 
could be significantly reduced by further trajectory shaping, or by intro- 
ducing constraints for this quantity. 

In Figure 29 is shown the trace of this representiVtive trajectory in* 
the Mach-altitude plane, at intervals of 8 seconds. 
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FIGURE 29. representative TRAJECTORY IN THE MACH-ALTITUDE PLANE 



Energy Variations During Zoom- Climb Maneuver 
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Other interesting results can be prepared from the output data 
tabulated in Table V. These relate to the kinetic and potential energies 
and their rates of change, and a computation of these variables has been 
carried out by W. Page, the NASA technical monitor of the project. The 
results obtained show how the energy components vary during the maneuver, 
and help to explain how energy is lost during the transfer from low to high 
altitudes. If the thrust could be varied during the zoom-climb so as to 
cancel the drag, of course, the theoretical maximum final altitude could be 
reached by converting all of the initial kinetic energy to potential energy. 
However, the thrust capability varies with Mach number and altitude in a 
different way than does the drag, which also depends strongly on angle of 
attack. In addition, the final velocity must be high enough to permit 
aerodynamic control of the aircraft. Consequently, of the total initial 
energy available, a portion remains as kinetic energy as the aircraft 
passes over the top. The velocities at this time vary from 550 fps to 
1150 fps, depending on the terminal dynamic pressure constraint imposed. 

The remaining energy loss is equivalent to only about 5000 feet altitude, 
and this is proportional to the time integral of the difference between 
thrust and drag during the maneuver. 

The energy rate has been derived in equation (200) as 


dE 

dt 





which is the sum of a "thrust energy" rate and a "drag energy" rate. In 
the energy state approximation, this is assumed to be equal to the sum of 
the potential energy rate and the kinetic energy rate 

^ = A + A A ^p , ^^k 

dt dt dt dt dt 

and these derivatives can be accurately approximated by finite differences 
read from the discrete printout in Table V. The result of the computation 
in a typical case is given in Figure 30, which shows in component form how 
the components of energy rate vary during the maneuver. Data for this example 
is given in Table V(e). 

In this example, the final dynamic pressure is 20 psf, and the curves 
show that the energy remains nearly constant during the maneuver. That is, 
the thrust and drag energy rates are nearly equal and opposite, as are the 
potential and kinetic energy rates. Of the total energy loss of 17940 ft, 
about 35% is due to the imbalance of thrust and drag, while 65% remains 
in the form of kinetic energy as the aircraft passes over the top. 

A second example of this type is shown in Figure 31, and in this 
example only the initial velocity differs from the values assumed in 
Figure 30. The Mach number is initially 2.1 in this case, for which numerical 
results are given in Table V(j). Here the total energy "loss" is 23430 ft, 
of which 52% is due to the imbalance of thrust and drag, and 48% remains in 
the form of final kinetic energy. 
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FIGURE 30. ZOOM-CLIMB TRAJECTORY ANALYSIS, ho = 50000 ft, Mq = 2.0, qfi^al = 20, 110% thrust 
Altitude reached 90034 

■i 

Initial Theoretical Altitude 107965 (at starting condition) 







The time -variation of the specific energy is shown in Figure; 32 for 
the two trajectories just discussed. It is seen that in both cases the 
energy increases slightly during the dive maneuver, when the thrust 
exceeds the drag, after which the energy reduces to its final value. 

These transients help to explain why the final altitude depends upon the 
initial Mach number, and why there is no unique answer to the question: 

What is the maximum altitude attainable for a given combination 
of aircraft, engine, wind and final dynamic pressure? 

The dependence of the final altitude on the initial aircraft velocity 
(or total energy) in the altitude maximization problem finds an analogy 
in the problem of minimizing the arc length from the point P to the 
line-segment boundary in the sketch below. ^ 


Nominal 

(example) 




Here it is seen that a numerical perturbation method may converge to either 
of the local minima (A or B) , depending principally on the nominal path 
assumed. Initial directions below the corner are likely to converge to B, 
and conversely for initial directions above the comer. By analogy, if 
the initial aircraft energy is less than the maximum attainable-, the process 
does not first optimize the "initial condition", but rather begins the 
dynamic portion of the maneuver, and this may lead to an inferior optimum, 
as described in the discussion under "PROGRAM VERIFICATION". 

When the initial velocity does not correspond to the maximum energy 
condition, as in the lower curve of Figure 32, the longitudinal acceleration 
due to the available thrust is extremely low. This means that a very long 
transient time (300 seconds or mote) would be required to increase this 
initial velocity to the maximum energy value shown as the upper curve of 
Figure 32. The numerical process in this case instead begins with a 
zoom-dive maneuver, and maximiim altitude is attained at about the same 
time as found for the upper curve, from the maximum energy initial condition. 
The final altitudes attained from these two initial velocities differ by 
less than. 300 ft and this is a measure of the "error" introduced by the 
existence of multiple extremals. 
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CONCLUSIONS 


It has been shown that a standard F4-C aircraft, zooming at safe 
operational weights, is capable of monitoring upper atmosphere pollutant 
levels. Generally, altitudes attained varied through the interval from 
85000 to 95000 ft, depending largely on terminal dynamic pressure, but 
the effect of improved thrust capability, of stratospheric winds and of 
reduced '..'ircraft mass on altitude performance have also been investigated. 

The peak normal accelerations on the optimal paths are high (up to 
4 g's), but it is likely that further trajectory shaping could reduce 
these peak accelerations without significant effect on the maximum altitude. 
Other problems which will require additional study include: 

i) Mechanization of the optimal paths for pilot guidance; 

ii) Analysis of F4-C handling qualities at low dynamic pressure; 

iij.) Shaping of re-entry flight paths; and 

iv) Effect of range safety and environmental acceptance of high- 
boom paths at prospective sampling sites. 
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TABLfi V, COMPUTER OUTPUT DATA 
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Run 
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Mq 
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a. 

1 

0 

0 

8.7 

2.00 

37500 

88191 

b. 

2 

0 

0 

20 



87167 

c. 

3 

0 

0 

40 



82780 

d. 

■ ■■ 

4 

+10% 

0 

10 



92645 

e. 

5 

+10% 

0 

20 



90033 

f. 

6 

+10% 

0 

40 



84928 

g- 

7 

+10% 

150 

10 



94625 

h. 

8 

+10% 

150 

20 



92612 

i. 

9 

+10% 
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. 40 

2.00 


87403 

j- 

10 

+10% 

0 

20 

2.10 


90371 

k. 

11 

+10% 

150 

20 

2.10 

37500 

93458 

1. 

12 

+10% 



0 

20 

2.10 

- 
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t.bfttft- 
> , ' 7 <J i ft 

I ^ < J » V 1 ft - 

ft 

. ^ 7^4 ft - 
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1 o3t .fthf* 
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ft?3 51 JSl -r.fr — 
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54 P 3 .ftoa 
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— ? T 5 e 7 - 5 ftO - 
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— 756 ?. 2 ftft - 
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-2557;a6«- 
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- 25577064 - 
2557, aou 
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- ?ft,2070V 
28.20706 
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-28771^79- 
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1 ,OQ9«4a 
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2 .‘ 24 eM 
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?,06'4259 
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2,094465 
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1,956638 
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1 ,ft032?4 
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l.7f«560 
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T , 069088 — 
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1 ,086689 
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1 , 0652P0 
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b'nsRoo? 

6'6aS60O 

T 661 TOl 
— 1R'pp? 673- 
7 0«iSnft1 
7 ?1 s?57- 
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3,9607 o5«E-02 
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1.95o20o9E-02 
■~3,95501o9fc-02- 
3,9598101E-02 
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4, l65lS36t-02 
— -4,28255 o2E-02 

4,O854484t»02 
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— 6,2159571E-02 
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— 6 , |05i 3Y6t-02. 

5, 43?94lOE-02 
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5,4090226t-02 
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4,59/6bA9E.-02 
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1 ,95/75o9E-02 
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I ,651601 U-02 
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-—319595b — 
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3174287 
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1172365 
- 3172550 


— CYNPP 

681,9943 

690,991a 

708,0266 

- 73 3,9062 
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1031 ,160 
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-15,21101 
12,33091 

-10,7l-.9o 
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10, 1479C 
10,14479 



ORIGINAL PA|GB IS 


Case Fj^*07 


*TnP III ATHOSPHtRIC TMAJECTORY optimization PKOCPaH VfcBSIOM J.OO JS MAY 1972 


STir.F 1 


ZOOM TO GIVEN 0, rORRECTFO OATa 


CYCLE 15 


PAUE 178 


VG 77 F — 

1<» JN.OOO 
ITaq.SPS 
1T'j5,6 05 

?r\^u,a7u 

- ?rU8,uib 
?ot‘', 7 «a 

,007- 
211?, ifeS 
-“21 51 , o«45- 
21 46,*! J9 

— 21 5*>,feSN 
215«*,U0» 
?lu7,hiS 
2122.552 

1 ,75« 
2021 .443 
•.“UP,!;)? 
1(»70,1=>9 
1 784 , S4? - 
I oOa , 971> 

- 1005, U4f- 
ISJ 1 ,?J0 

' 1-.2),0”S 
1311 . S4S 
!?3'>, 2? 
nS3. 548 

— I 0 7S, 4 I 8 
1007,258 
9S0,844o 

*>0/, 781 0 

♦'*0,2521 

86'7.2S71 

— 8^9^287 A— 
8 e»'>, 2 S 7 « 
8fc9, 1034 


MnCTF 

50000.00 

40017,92 

U78m7,5*7 

«o?uo,e5 

4“m48,*>0 
~B7*7a,5b - 
48952,32 
458 30,20 
4as<Jo,?5 
«. 5?So, 3(i 
« 18 S?, 43 

-0004?, Oti — 

39092 , .i 7 
37892, O'i' — 
37001 ,05 
“ 38870 , 75 1 — 

37178.71 
-3MM9,83 

41115,10 
-«45e8,8« - 
a*»*,8o,87 
-55181.71 — 
57858,57 
-62“5u , 07 — 
88808,54 
71880,37 

TiiOpO , 1 8 

78585.52 
« 1 589,50 
* 597 a,o 9 
8bl09.72 

87787^81 

89001 ,Po 

89751 . 1 I 

90031 .52 
-900 31 .52 — 

90051 .52 

90055.71 


CAM70 

0 , 

-1,21 7508 
-2 , 482504 
-3,715211 
-4.951213 
-6,147703^ 
-7, 289025 
• 8,2809>J2 
•9,017470 
-9,485155- 
-9,829102 
-»9, 407230- 
-0. 075123 
—•f, 18082^ 
-4,59595a 

, 5980510 

6,482013 
--1 3, 81897- 
21 ,86448 
29,452^8 
35,55163 

- - 40,26 750 - 

45,40322 

4S, 18882- 

45,62147 

44,97872 

-5,57730 

40,8u^23 

37, 48207 

5 5. 29 7 55- 

28,25540 

2?, 381 55 

15,70001 

6,415490 

,7465045 

— , 7485843- 

, 7465843 

0 . 


— V1T7P 

3485,825 

3473,954- 

3484,414 

- 549 7,016 
3511.729 

“3528.424- 
3548,686 
5568,644 
5567,091 
3808,095- 
3623,878 
-3o48, 581- 
3685,2-4 
3676 , 360- 
3684 , 452 
-3878,287- 
3645.372 
35*5,687- 
3*.67, 125 
-5565.981 
5239,088 
-3115,095- 
2990,574 
•2895,803- 
2804, I 77 
2727. 155 
2881 . I -5 
2599,035 
2545 , 79o 

- 2501 ,834- 
2486,029 
2436,514- 
2416,888 
2408,751- 
2401 , 755 

■240 1 , 755-- 
?40l , 755 

2401 ,848 


9CT7M 

0 . 

—1,275942 

2,555850 

3,640690 

5,135511 

— 6,459281- 
7, 755127 
9,077694 
10,41587 

— 1 1 , 78246- 
15,12411 

14 , 4 V9yftr 

15,66715 

— 1 7,2m750 
18,69618 
20, 1 1245 
21 ,51288 

-22,87398- 
24,15790 
25, 33168 
28 , 3«685 

— 27,55“17 
2«, 13328 

— 28,95053- 
29,87294 
30,54655 
30,99455 
31,81915 
32,22851 

— 52,62064- 
33,40ol0 

— 33,9m49« 
34,55944 
15, 1 511 S 
55.70U2 

— 55^70142- 
35,70142 
35 , 7584(3 


— a-aCh — 
1 ,999644 
-2,008853 

2,020095 
-2,034198- 
2,050901 
2,C700 50“ 
2,091255 
2,113901 
2,138985 
-2,159982 
2, 182025 
"2,2ui 73«— 

2,21 7426 
2,2277/8 - 
2,229588 
2,2 10458— 
2, 192320 
-2,150383- 

2,088104 

2,013564 

1,95U72 

l,*-339i 

1 , 75Ce70 
-4 ,65831 7- 
l , 361065 
1 ,40o9 1 H 
I ,37'>2^9 
1 , 2 8 0 S 

1.191584 
-I ,107109- 
1,032549 
,9710001 
,9249164 
,8954649 

,6658176 
,88361 76- 
,6656176 

a 68 56568 


- AM456 1 

17502.76 
lTfcC6,45 
17674,90 

n e 1 c 7 , 4 7 

18709.76 

-1946a, 2^ 

20“24, 79 
21493,25 

22761 ,55 
244C2 , 4 7 
26066,82 
-27 7 79,14- 
295r /,5u 
3!il.*5,6<.- 
325?-, 53 

- 3 3 C' 7 c , 2 5 
32342, 75 

-301 72,87 - 
267r6, I 7 
2232^,-3 
I6lt3,96 
14332,83 
11064,50 
-8392,595- 

6149,985 

452o , 7»-3 
5364 , 2‘- 1 

0. 

0 . 


PAGE J8 


TABLE V(e) Run 5 



CtSE AjM40? 


CYCLE IS 


?T*GE - 1 


ZOOM TO ClV^N 0, rOouECTCO PAT* 


AL«’XP 


E77P 


fr, ian/JP7t*0? 
S,o?l lS7«t-02 
S,PP*»«Pr'«t-0? 
S, C^‘^e5a U -> 0 ^ 

5,01«V«b5i-o2 

a, 0a«i4i 1 7k. 0^ 
5,0«^e«'S‘>t.o^ 
S, 1 7oS7S0t-p2 
S.n‘>SSSlfc-02 
S, J70‘^Sl''»t««)Z 
5, 7<f5JOOdE-o2 

b, 70?ZOOtk-oa 

«,O^67'»b0t-0^ 
— ,iiab7os 


?,^??'‘37 

a. is'7?HS 
?.0‘>l 310 
Z.ftIPP®® 
1 ,V^*)'»V7 
J ,«O7?07 
t-BTTOo? 
1 .'>^^S7 J 
-t,0'«'iOi)7 

1 ,V5<J0U5 

2. I ai«b7 
^,uuS'J^7 

« , 0S7 ?60 

6,0?5^S2 
6, I 7«9S0 
6,0??X7 J 

b, 1T90Z6 
■4>^q7.C,4^ 
S,bSb07 % 
S.?'Oi 

« , 79on I b 
-4»,<i?bS!»» 
0. 1?7i«7 

J.oSat IS 
J.S9J0J\ 
3 ,SSP0«7 
3 , S » 

A , 3 

.3,«iAS«;b3 

S.USS'-bJ 

3,4f 


IbfiOP.Ub 

171 5?, 6b 
tTSO?.*37 
JAi?ft,7? 


|o7 70,S'» 
ZO^O?,*!’ 
ZZOZo.lZ 

Z3337,77 

?si Sb.S'J 

?PU 03 ,Sb 

3pooa,^^ 

3l?3S,6« 

3l7rtH,br 

3l?S?,6S 

2bl hS,u2 

?1«P3,7S 

17«S7,‘>2 

1 3«b9,30 
10300,75 
- 73 PS, 70<5- 
S2M ,0?l 
37P9, 37^ 
2770, 2bl 


lbSS7S9 
1 9 u09O4 
21 13309 
a2b**70b 
2319015 


?3«2bl2 
?31 9704 
?2'>4 ?o4 
?^9SC-<3 

237Pb90 

?4tt40b> 

2933212 
2462123 
20 7bo 34 
24201 62 
2572623 
2372H23 
2372P23 
23700 35 




rKTXtrV'ro'CTj-.A'i V^RSIO** 3. UU pay 1972 


n 



PASS 2 


pace 179 


— CO »v?B 7c tOPt7^ OTKPP 

4,03076296-02 -,bS9-22i 3973ob7 b6l,99«5 

■ 4, 02ol 16bE-02 — -, 6990227 39tf7c00 6*'C, 722S 

i,0069u62E-02 -,bul26Sa 3S01230 707,0302 

l,'’97i2P-*k- -2 -,63<>f706 3SJ9SS1 75l,S9cC 

.l,VHu2b62E -02 -,6405046 352/492 7bS,22S2 

^l,9 7 396 3^t-02—-, 65229 TH 35 39n3C 606,910 3 

i,9^4^83iE-|1^ -,6710206 3551 162 663,6360 

l,96iPH75k -0 2- -, 7l 31553 35»07lp 930.29 7 6 

;.,9669696E-P2 -, 7939767 3567 1 3a 1006,775 

-:.,97j2379E-02 - -,6 79e09} 55 70727 IC’d.beO 

] ,9571/1 lC-02 -,9S049 o7 3571054 1199,005 

,966 7o39fc«02 — -t • 1*15063 35oo572 1 306,01 7 

A ,003373bt-02 -1,25 j9SQ 3555965 1915,192 

i^05rtl006t-02 - -1,551940 3S392<ii l5lO,600 

A,l596593E-02 -2,025967 3519582 15/9,228 

•97379 1 3 70E-02 --2,71 52 5 7 3««089S 1526,9 39 

«,«32565Pf-02 -2,729051 i«93C37 1513,990 

-9,899258«*k-82 — -3,22 722* 34 0 3685 1 359, a5« 

•.t,25l«257E-02 -5,19C22* 3359e02 1 1 37,465 

V ,322000 ?E -02 --2,6 33793 3 326427 *96,o5*>l 

6,56791. ;St-02 -2,1691S2 33064^8 6/7,6835 

5, 76;6M60e-02 • -1 , 706 7 30 3 29i»»8- 9 97,6 646 

3,75779?lt-02 -1,262100 3285417 S59,J6b6 

-5,6 7 36 7596-0 2— ,914660? 5261 |9« 25 7,9 3fcS 

5,550901 7E-02 -,6558228 3278636 185,3149 

5,39 799|>E-f,2 -,4,9i2gs 327 Tw5S 159, ''738 — 

5, 1975396E-02 -,3592788 3277079 96,10811 

9. 9651 125E-02 -,2450153 3268464 72,12047 

U,8Ul6OOnE-02 -,169<1 p 5J 32oll7l 59,09631 

-9,60 319986-02 — I 48 796 7 3256 1 27 9 | ,2 7059 

9,l69S962f -02 -,116298a 3252650 32,94900 

3, 4 323029k -1*2 -9,6«26 3!8t-02 3250290 26,54; 28 

2,487C6p4f-i? -8,26993501 .02 3240.-58 22.75458 

2,1652662E-02 -7, 3 315257k-P2 329/905 2P,«.i)692 

2,064623/t-C2 -6,»8dl475r-P2 324716a l9, 81510 

2 ,064623/f “02 - - 6,8\.*6|975k -0 2- 324 7 1 6. 1 9 , 8 l 5 I 0 

2 ,069623/1-02 -6,8t«19rS£-02 529 /160 19,61510 

2,0t23908k-02 -6,879i595f-P2 324/099 19,80566 



I 

i 


*TnP tit *THt)SPHti»ic trajectory OPT I**1 2*T PROGRAh vtRSIO^ 3,00 15 h4v 

CaSP raWAO? Stage 1 CYCLE 15 PASS 2 PA« 175 

Zoom to given o, corrected data 


TI*^E 

1 0 . 

3 noof'flo 

5 Ib.OOftOi- 

^O.OOOOfl — 

y 2«.0florto 

0 ?'4,00'>nft 

R 3?,nnono 

10 ^''.OrtftOO 

11 ‘*0.00000 

12 40^00000- 

13 4?.0O0'‘O 

1“ se.ftftooo 

l5 •*(«.•» 0000 

l» »'u,0''0ft0 

l7^ ♦'O^oOftOO 

I** Y2^3/>nno 

20 70.00000 

?I '•O.OAOOO 

22 *‘-•,00000 

21 *‘'1.00000 

2** • e , OOOOO — 

‘>*'.00000 
2 P 1 ■ 1 . 0 0 n 0 4» 

2 7 1 0 , , 0 o 

f 1 f , looo 

2R 1 H , 1000 

30 1 1 o , OPAO 

il ' I 0.72A0 

32 110,7^70 

33 110,9^70 

30 110,7570 

15 110,0570 


VOT/P 

19Ja,O00 
lOUO.OUO 
I 9 7 1 , K 'i 9 
j qgA 

20?4,7*.h 
- ?OS0,?7O - 
2060, ?hS 
2072,727 
20*-R , 760 
2oS 5,R«2 
2020, ?,52 
- JOQ^ ,09u- 

1 9a2 , *R*> 
1*90,204 
1«3?,021 
1 7^7, 5U2 
1 700,05 J 
1651 ,750 
1562, u**! 

1495,691 

1 426,592 
1 i 6 I , 4 U 0 
1 ?05,6H0 
-1235,650- 
11*2,345 
1 140,6^7 
! 1 (.9 , 7r>0 
I 0 9 r. , A 0 9 
1 090, h09 
1 0<»4 , 4 J9 
1064,439 
1064,450 
1064,450 
1 0*4 , 4 50 
1064,450 


— **GC7P GAH7i) 91776 RG77*« A»**Ch 

50000,00 0, 3465,625 0, 1,999644 

-497 51.51 -S, 67151 4 3477,569 1, 274*06 2,<M424j 

469hs,65 -7,257977 3494,290 2,555625 2,0566t« 

-47777,23 •9, 716271 1515,697 5,644479 2,064595 

46559.55 -10,55599 3541,295 5, 145146 2,095^52 

-«a69n,s6 -9,t,2iMoo 3567, 1 79 6, 405465 2, 1 1 7*“'2- 

45665,27 -7,540065 55*8,240 7,000«50 2,154425 

42«51,u 5 *4,0 59954 3599,645 9,iS3955 2,141060 

“2522,09 -,2655176 3598,854 10.51591 2,158018 

42760,47 4,245462 3582,574 11 ,86887 2,121679- 

45761.45 9,36817a 3547.172 13,20142 2. 0965*4 

-454 28,29 14,6 8888 34 92,001 14,49 5 74 2,0So7 55 

47764,39 19,65614 3422.130 . 15.72730 2,046970 

5(159, ,62 — 25,49640 5549,037 16,89741 1,952557 

55749,95 26,62479 3272,521 18,00432 1,692435 

57124.55 29,19996 3|9i,569 19,04835 1,82562“ 

60601.57 50,75578 3115,119 20.05297 l,75ell6 

64 04 2,8 4 51 ,5 3952 50 44,989 20,969 4 6 I ,6*5560 

67367,06 31.25866 2979,007 21.86«56 l.6l«o07 

76S12.04 50,58 4 65 29 I 7.6 32 22,724 1* 1 ,5-86 45 

75455.42 29,52194 2061,790 25,55358 1,473637 

76080,50 27,44155 28U.662 -^24,35774 1,406295- 

7»o 5«.27 20,98471 2758, 97o 25,15817 1,556550 

-80 45 5.35 21 ,89 640 2715, 752 25,89 755 1,2/4 532- 

82112,98 18,19655 2680,751 26,64089 1,2/1553 

83402.54 15,94917 2655.576 — 27,37182 1,175682 

80508,23 9.195674 2655.611 28.09572 l,19lc«2 

8uai6,67 4,007688 - 2621,215 - 28,80925 1,12079.1 

*4816,67 4,007668 26?l,2l5 28,60925 1,120790 

64 Q 2 A,aA 1 ,245752Cfc-02 - 2»»l6,6l 1 29,35050- 1,114189- 

*4928,86 1 ,2457520t-02 2616,611 29,33050 1,114189 

*4928,88 4,02351551-06 2616,602 29,332H 1,114180 

«4'’28 .*8 4,02351 55t. 06 2016,602 29,3521 1 1,1 19180 

8<j 92A,8H 4,02551 55k-o6 2616,602 29,5521 1 l.lUl*© 

*4928.86 0, 2616.602 29,35211 1.119180 


-ANAS* 1 

17502,76 
17596,07 
18567,69 
19599,26 - 

20975,47 

2 / 3 6 7 , 6 V 

2.5660,04 
24851 ,92 

25207.75 
24862,63 
23654,91 
21657,42 

19J27,2^ 

16717,56 

14209,15 

11910.76 
9822,589 

7961. 7p2 

63j6, 772 
5025,496 
4050 , 9 S 4 
0 , 

0 . 

0 , 

0 , 

0 . 

0. 

• 

C. 

0! 

3t 

0 , 

- 0| 

0 . 



TABU; y(f) Ru n 6 


*■ I 
» • 


»TrP Til ATM0SPHE3IC TWaJECToO'' UPTIHIZ4T IC»< ■'Knc«i“ VT»5I0Ni 3,00 IS “iv 1Q72 


CiSF r<jM*oT 


STACt t 


ZOOM TO RIVE'' 0, cPPSECTfO 0 *Ta 


CVCLE 15 


P*UE 176 


AVAS 5 
1166,60') 
I 16'J, 

I IhS.SRB 

I t 60. 

1 tSP. JRO 
1 1S7.7SS 
1 1 Sh. 670 
I T\7 

1 16 <». 60 T 
I ISO . O ?! 
1 t 3<J, TS6 
I l<t7,'»T«> 
1 , hOS 
1 ^ 

1 I '< , 7 ? 7 
1 (1 \,970 
1 0 1. T65 
I <i^.«73 
1 ■4^,063 
1 <JP. 1 71 
» II .OS 7 
t -11 .PS 7 
1 '!t .PS7 

» 41 '.os 7 

I'll , PS 7 
* Jl .PS7 
1 '1 .OS7 
1 J 1 .PS 7 
1 'I 1 .PS 7 
I 11 , OS 7 
1"1 .OS7 
101 .OS7 
1 J I .PS7 
1 '1 1 PS 7 


At PHft 

»,3P*S1 S3E-03 
,T«176PP 
,6U4S?7<1 
l.SS0?5P - - 
?.S»1 06*1 

-3.a30*;iP 

5.753050 

а . 20 nP 6 P 

u . 06S« 1 a 

5.000566 

5.0051 Pi 

-5.030571 

S.70P37S 

5 . rsoSPi — 

б . 1 ? Po,»P 

6.61 <717 — 

6,0?01 ns 

-6,'16070S 

6. J31 <26 
6 . f * PS 5 Hl -- 
6.675605 

7.01 67^6 

7 , ;> 0 P 1 '4 6 

t ./ s '- S ?? 

6.5PS060 

6.;»0S772 

5. 757P67 
U ^ gono^S 
0 *PP5P;>S 

o ] 6 S 1630 

0.651 6 <0 
O.6S0S60 
0 . 65 (I S 6 0 
0,650560 
O . 6S0S6O 

» t > o 


— rrnp — 
1 6*00, «» 
17066, oP- 
1 7505,00 
-1 5050,0(H 
?035S,OP 
-?i 701 ,3r- 
?51 1 1 ,?*> 
20077, to 
20006, 60 
20162. n^- 

21055,61 
"?1 1?5.hf- 
1 5570,65 
16200, PI 
13667.25 
11205.10 
25S6’.650 
-«>oi 0,515- 

5 U 00 . H 66 
<»21«.ol5 
1167'.651 
— 0 . 

0 . 

— 0 .- — 


.1 ,07lS255h-02 
•1,52561 7ot-o5 
7,1160521 lfc-01 
1,56310 J2fc-02 

7,25P0506fc-o2 

— , 1 000661 

, 1102076 
- ,1206215- 
, 120352O 

—,1506041 

,1512165 

— 151PP75 

,1 755751 

,1561P52 

,2057«67 

,2115121 

,21 74162 

--, 2004PP7 

,2511000 

,2021111 

,2P51hP7 
, 1121060 
, 1626242 

- , 1602061 

,1721560 
, 167U2P7 
, 1S24068 
,1122277 
,1122277 
,^021575 
,20215/5 
,2021220 
,2021220 
,2021220 
,2021225 


— t ; i > 

3,8662505F-02 

l,«llou0SE-02 

1 ,8075'471E-02 
1.01055111-02 
4,05560071-02 
-«l,21S756oE-52 
0,1106501E-02 

- 4,«5470?»t-02 

а, p413527E-02 
«, 750401 lk-02 
O,7701116t-02 

■5,2056544f-02 

5,0902-»lPfc-02 
5,161 700PF-02 
5,01b460lE-02 
•5,«421 121f -02 
5,505?560F-n2 
-5,642736lF-02 
5,022lJ06t-02 

б , 5?45lOot-02 
6 .41071O0F-02 
6 ,P 1 60«6l^ -02 
7,275l6U5t-02 
6 ,P1 Mi551fc-02 

6 ,60u5'474t -02 
o ,54P10O6t-02 
6,10 160 /of -o2 
5 ,44 lSS20f -02 
5,40l512Of-02 
5,lS64ll6f-02 
S,l55'Jli6t-02 
5 , l5/66l2f -02 
5,l5/66l2f-o2 
S,l5/66l2f -02 
5,l5/h6l0f-02 


5N75/C 

, 1055104 
-2,1 15«807F-02 

-, I 2250UQ 
-,4022376 - — 
-,0700251 

- -1 ,46ttSl 1 

-1,715650 
- 2,050445 

-2,008764 

-2,44u7Ht, 

-2,113604 

- -2,U66lC0 

-2, 042506 
-1,761672 - — 
-1,571607 
-1.10051 7 

-1 ,010466 

-,004 104 4 

-, 7110554 

-, 6610662 - 

-,5164155 

-,47240S6 

-,416/056 

-,l4lll4S 

- , 20006*2 

-.2565621 - 

-.2221215 

-,1551525 

-,16Si525 

-, 1 701546 

-,1701546 
-,i7omiQ - 
-.1701410 
-,1701410 
-, 1 701434 


E 5 PEF 

3*71667 

- - 3402106 

1510352 

1526200- 

1536513 

1535650- 

1532567 
- 1510116- 
1501114 

— 1451203- 
I4o061 0 

1436221- 

3415550 

3404800- 

1507124 

1 - 

1591 O 20 

-1177445- 

1172370 

-- 1166940 

1161652 

- 115-015 
lllO'.So 

1 127152- 

111 7011 
1110125 
51C-21 1 
1200215- 
1200215 

120605} - 

1206051 

— 1206044 
1206U44 
120-044 
12O604U 


— C 7 S 4 P 

681,0041 
700,8002 -- 
741,2208 
-«65,42l6 - 
800, Util 
-07e, 7002 — 
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-.691 1739 

34,'»<, 1/ 

205,-9/9 


t 

1 

2(> 

1127.616 

5.255971 

407 1, |Oft 

,2166«55 

5.2700161F 

-02 

- , 485*4 1 1 

3 426066 

1 48 ,2» 1 4 




27 

1 127,36ft 

0.804617 

Of 

,2127147 

5,091 *1201 

-02 

-. l40*4ft5 

34229*1 

107,8790 




28 

1 1 27,36ft 

0,418051 

0. 

,2112782 

«,919652if 

-02 

-,?47*5i 3 

3 ft i (• 4 a 5 

7*.9«lo7 


f 

■ 

2’ 

1 1 27, 36ft 

0. 1 2011 0 

0. 

,2129122 

1 , 7*02 7'46F 

-02 

-.1 *60132 

3401999 

56.9984* 




30 

1127.36ft 

1.6601 76 

0. 

,2221 177 

4,6*697 763. 

-02 

- , 1 4 7o io2 

3 1966*4 

40,76173 




31 

1127,164 

3,661751 

0, 

,2302771 

4 ,49866427 

-02 

-.1192201 

1 19187© 

35,0*594 


0 — . _ 


32 

1127.36ft 

1 . 0 0 4 7 1 5 

0. 

,23771 12 

4.24519*37 

-02 

-9,986567;>f-02 

31t*r|2 

2 * , 4 8 ; 1 0 



|r 

33 

1127, l6ft 

3.560517 

9. 

,2422832 

3.91250 7S3 

-02 

-*,61 '131273 -02 

3 366S6© 

2«. 1 173* 




3« 

1127, I6ft 

3. 5201 12 

0, 

,2452510 

i,2128j2li3 

-Oc 

-7 , 7 39,i«fc7t-i»2 

31*4921 

21, -6522 


/ 


35 

1 1 27.16ft 

3^ aoftooO 

0. 

,2451244 

2, 77*29723 

-92 

-7,2'>.’2764I -02 

11M7S6 

20, 14824 



1 

36 

1177, l6ft 

i^oftaoo 

0. 

,24512114 

2, 7 7«?972» 

-*2 

-7,25227843-02 

118*75© 

20, 1 4624 



1 

37 

1 127,16ft 

1. 084 1 VI 

Dt 

,2440617 

2,71 305951 

-02 

-7, 1427805E-02 

336M93 

19,931 i9 


t 

1 


3* 

1127,368 

3.4*4191 

0. 

,2440637 

2.7110S95F-02 -7, 1 u27*05f-02 

3181193 

»9, 91119 

3’ 

1 127,3«ft 

J',4*4191 

0, 

,2440637 

2,71 30595F-O2 -7 , 1 42 78o5E -02 

316U91 

l«,911 19 

oO 

1127,368 

3,4*4180 

0. 

,2440629 

2,71 10458F-02 -7, l42752*k-02 

3161193 

»«,911 1 7 











f 

t-- 

f 

iu 

/ 



CiSE r<t«A07 


STAGE 1 

C7CLE is 


PASS 2 

PAGE 175 



2POM 

TO given U, 

cn«9tcTE0 data 







T?mf 

WG777 

HGC7f 

0AM70 

VI77E 

RG77N 

amacm 

amas^ 1 

1 

0. 

2011,000 

SOnOO.OO 

0. 

5500,625 

0, 

1,999800 

l7iC2,76 

2 

«, 009000 

202‘>,0i'» 

«960S , \h 

•9, 506667 

5552,177 

1.325039 

2,012157 

17629,56 

i 

• , "ooooo 

20ii«, 1 'iO 

c»0ia,12 

.7, 009085 

5569,592 

2,652682 

2,029087 

18077, 66 

a 

t?,00000 

207»-, iu« 

0 7U69,<1 t 

-10,69907 

5592,259 

3,989161 

2,051 559 

19770,00 

5 
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1,951012 

15621,69 


1 A K) 

AA‘^ OCOOO 

1«57,S90 
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26,01001 
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29, *.8197 
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2« 
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78501,59 

27,59850 

2860,906 
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0. 


25 
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1507,772 
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20,88079 

2812.500 
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0, 


26 
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29,71099 
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in 
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1155,099 
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0. 
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2687,052 
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0. 
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2687,006 
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0, 
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0. 
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A7ao A ,SS 
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2 6 8 7 , 0 9 6 
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0. 

5 * 
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26*7,006 
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0, 

C 
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0. 
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a 
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5 
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1001,592 

T 
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a 
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• 1 ,829190 
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10 
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11 
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• 2,'.'10«29 
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1« 
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,1606360 
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,20608'»8 

5, 5628751E-02 
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,26 35868 
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0. 

,3u3l030 

6.9'l%82?2l - 02 

— #91^6999 
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0, 
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0, 
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3953921 
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APPENDIX A 


PAST COMPARISONS BETWEEN PREDICTED OPTIMAL PATHS 
AND actual flights PATHS FLOWS 


The success of the variational steepest-descent method in solution of 
aircraft performance optimization problems is evident from the strong 
support given to this technique by a series of contracts let by leading 
Government research centers concerned with this area. The reason for this 
support is clear when performance gains obtained are examined. Figure 1 
presents the 1962 time-to-climb record flights of the McDonnell F-4B aircraft. 
Figure 2 illustrates how closely these paths follow the minimum time ascent 
paths predicted by the References 1 and 2 program. Figure 3 provides a 
comparison between flight handbook performance estimates, a minimum time 
climb obtained by the References 7 and 8 program, and an attempt by Marine 
Col. Yunck to fly the predicted optimum. 

The predicted optimal path and the path flown by Col. Yunck both produce 
a 23 percent improvement in aircraft performance over the flight handbook. 
During the Cuban crisis of the early sixties, results of this type were 
produced routinely from the References 1 and 2 program to aid in Air Force 
readiness studies. It should be noted that unlike optimization studies in 
other technology fields, these performance gains are obtained without vehicle 
modification. To obtain these performance gains while retaining flight 
handbook methods would have required a 23 per cent increase in aircraft design 
capability, several years’ effort and several billion dollars to replace 
an existing fleet of aircraft which could achieve this capability simply by 
being flown in the optimum manner. This one example serves as a lasting 
example of: 

1) The high cost associated with an oversimplified approach to 
performance optimization; and 

2) ' The insignificant computational cost of adequate performance 

optimization studies for production aircraft when compared to 
the resulting payoff. 

Further details of these F-4B performance optimization studies may be 
obtained from Reference 5. 
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